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INTRODUCTION 


Our  primary  objective  in  this  research  program  was  to  use  experimental 
atomic  beam  scattering  techniques  on  carefully  selected  systems  to  increase 
our  understanding  of  the  detailed  collisional  processes  occurring  in  three- 
body  (atomic  plus  diatomic)  systems  Involving  multiple  electronic  states  with 
many  competing  product  channels.  We  have  used  a  sophisticated  atomic  beam 
scattering  apparatus ^  to  study  and  characterize  various  important  reactive  and 
inelastic  channels  in  collision  systems  that  are  complex  enough  to  challenge 
the  most  advanced  experimental  and  theoretical  techniques  available.  By 
systematically  varying  the  projectile,  the  target  molecule,  the  collision 
energy,  and  the  process  investigated,  we  have  been  exploring  a  wide  variety  of 
interactions,  searching  for  insight  into  the  competition  between  various 
reaction  pathways  in  relatively  complex  systems. 

Earlier  work  on  two  different  classes  of  problems  had  laid  a  foundation 
on  which  to  build  our  efforts.  First,  numerous  molecular  beam  experiments*  at 
energies  near  thermal,  accompanied  by  theoretical  studies  ranging  from  clas¬ 
sical  trajectory  calculations  to  quantum  dose-coupling  techniques,  had  pro¬ 
vided  a  great  deal  of  insight  into  both  reactive  and  nonreactlve  collisions  of 
an  atom  with  a  diatomic  molecule  in  the  special  (and  important)  case  where  the 
predominant  interactions  were  confined  to  a  single  electronic  potential  energy 
surface.  Second,  a  large  body  of  knowledge^  had  been  acquired,  both 
experimentally  and  theoretically,  regarding  electronic  excitation  processes, 
energy  transfer,  and  ionization  in  colliding  systems  composed  of  two  atoms  (or 
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an  atonic  ion  and  an  atom)  where  two  or  more  electronic  states  were  involved. 
Systems  that  conbined  these  two  features  (three  or  more  atoms  and  multiple, 
competing  product  electronic  states)  had  received  little  attention. 

Although  many  complex  trlatomlc  systems  are  amenable  to  detailed  investi¬ 
gation,  most  would  yield  results  so  complicated  as  to  be  completely  intract¬ 
able.  Only  a  carefully  selected  set  of  reactions  could  be  expected  to  yield 
any  true  insight  into  the  processes  occurring  on  multiple  and  competing 

i 

potential  energy  surfaces.  When  our  work  began,  the  greatest  success*  In 
characterizing  a  complex  trlatomlc  system  had  been  achieved  for  the  interac¬ 
tion  of  fast  alkali  atoms  (M)  with  halogen  molecules  (X2).  A  detailed 
examination  of  ion-pair  formation  in  this  group  of  reactions  had  produced  a 
generally  clear  picture  of  the  dynamic  properties  Important  for  describing  the 
Interaction  on  the  two  relevant  intersecting  electronic  surfaces. 

We  Intended  to  build  on  this  work  by  choosing  systems  (rare  gas  metast¬ 
able  atoms,  Rg*.  colliding  with  halogen  molecule  targets)  that  had  analogous 
ion-pair  formation  channels  as  well  as  a  competing  set  of  processes.  One 
major  competing  process  in  these  systems.  Penning  ionization.  Involved  an 
ionization  continuum;  such  continua  had  not  been  studied  in  any  great  detail, 
and  they  were  a  challenging  research  frontier  both  experimentally  and 
theoretically. 

Initially,  we  investigated5  ion-pair  formation  in  Ar*  +  I2  collisions  and 
searched  for  similarities  and  differences  in  ion-pair  formation  relative  to 
analogous  M  +  X2  systems.  The  similarities  were  striking  and  suggested  a  very 
strong  analogy  in  ion-pair  formation  for  M  +  X2  and  Rg*  +  X2  systems.  The 
differences®  were  minor  and  yielded  little  insight  into  the  interesting 
continuum  channel. 
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The  direction  of  our  work  then  changed,  as  we  began  to  examine  systems 


(Rg*  +  02  In  particular)  and  product  channels  where  competition  from  continuum 
processes  was  expected  to  be  more  easily  observed.  We  worked  closely  with 
Professor  Joop  Los  and  Dr.  Aart  Kleyn  (at  the  FCM  Institute  in  Amsterdam),  who 
were  simultaneously  investigating  M  +  02  collisions.  In  these  systems,  we 
observed^  some  fascinating  consequences  of  the  rapid  (subpicosecond)  vibra¬ 
tional  motion  of  the  traslent  0~  ion  formed  at  the  instant  of  ion-pair 
formation.  Again,  the  Rg*  and  M  projectiles  gave  very  similar  ion-pair 
formation  and  inelastic  scattering  results.  There  were  no  obvious  modifi¬ 
cations  of  these  channels  in  Rg*  +  02  caused  by  the  competing  continuum 
processes. 

It  became  apparent  that  the  continuum  channel,  even  when  it  was  the 
dominant  Inelastic  process,  was  not  easy  to  observe  Indirectly  through  its 
effects  on  the  product  distributions  for  other  channels. 

The  last  major  goal  of  our  work,  therefore,  was  to  develop  a  method  for 
examining  the  continuum  processes  directly.  To  do  this,  we  developed  and 
tested  a  delayed  coincidence  technique  that  allows  for  the  separation  of 
Penning  processes  from  other  Inelastic  and  elastic  scattering  events.  First 
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results  for  He  (2  S)  +  02  are  presented  here. 

The  collision  systems  we  have  studied  are  typical  of  collislonal 
interactions  important  in  a  wide  area  of  gas  phase  physics  (e.g.,  discharges, 
lasers,  excited  atmospheres)  that  require  Improved  characterization  and 
understanding.  As  more  detailed  reaction  rate  information  is  generated  and 
understood  for  a  few  relatively  simple  and  carefully  chosen  model  systems, 
comparison  and  extrapolation  should  yield  increasing  Insight  into  the 
mechanisms  important  in  the  most  general  cases. 


3 


i ua:« 


RESULTS 


Ar*  +  I? 

At  collision  energies  above  several  eV,  Che  inCeracCion  of  a  rare  gas 
mecasCable  (Rg*)  wich  halogen  CargeC  molecules  (X2)  has  an  imporCanC  reacclon 
channel,  ion-pair  formacion, 

Rg*  +  X2  ♦  Rg+  +  Xj  (or  X  +  X*).  (1) 

This  channel  resulcs  from  an  elecCron  jump  occurring  aC  Che  surface  crossing 
between  che  relaclvely  flac  Incoming  covalenc  pocencial  surface  and  an 
(asympcoclcally  higher)  aCCracdve  coulomblc  Rg+  +  X2  surface  ChaC  cues 
chrough  che  ini rial  covalenc  scace  aC  lnCermedlaCe  or  large  incernuclear 
separaclon.  A  similar  channel,  Involving  analogous  surfaces,  is  ImporCanC  in 
collisions  involving  alkali  atoms  M  wich  che  same  CargeCs, 

M  +  X2  ♦  M+  +  X2.  (2) 

For  alkali  projecCiles,  this  surface  crossing  dominates  Che  reaction  dynamics 

and  initiates  Che  prototype  harpooning  chemical  reaction  long  studied  at 
a 

thermal  energies 


M  +  X2  ♦  M+  +  X2  ♦  M+X"  +  X. 


(3) 
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The  second  seep  rapidly  follows  the  initial  electron  jump  because  the  strong 
ionic  M+X-  bond  can  form  at  the  expense  of  dissociating  the  much  weaker 
transiently  formed  Xj£  bond.  Even  at  higher  energies,14  when  the  projectile's 
velocity  precludes  new  bond  formation,  the  single-surface  crossing  continues 
to  dominate  the  dynamics,  and  the  major  Inelastic  channel  evolves  to  ion-pair 
formation. 

Similarities  in  electronic  configuration,  ionization  potentials,  and 

it  O 

reaction  properties  of  M  and  Rg  have  often  been  noted.  The  alkali  analogy 
is  a  useful  one,  at  least  as  a  first  approximation,  for  understanding  both  the 
ion-pair  formation  [reaction  (1)]  and  certain  collislonal  excitation  channels 
for  Rg*  +  X2  collisions.  This  analogy  is  an  Important  starting  point  because 
the  M  +  X2  reactions  have  been  studied  experimentally  in  great  detail  and 
because  the  existing  theoretical  models  for  alkali  reactions  can  be  used  to 
predict  the  corresponding  Rg*  interactions. 

However,  the  differences  between  the  two  systems  are  also  fascinating. 
Most  significant  is  the  large  Rg*  excitation  energy.  If  this  energy  is  larger 

than  the  ionization  energy  of  X2,  the  Penning  ionization  reaction 

Rg*  +  X2  -  Rg  +  X2+  +  e-  (4) 

can  occur.  In  this  case  the  original  covalent  Rg*  +  X2  potential  surface  and 

the  interacting  Rg+  +  X2~  surface  are  imbedded  in  an  ionization  continuum  that 

can  deplete  intensity  from  the  two  alkali-like  surfaces  during  the  collision. 
This  extra  channel  and  the  manner  in  which  it  modifies  and  competes  with  ion- 
pair  formation  add  an  extra  dimension  to  the  Rg*  +  X2  systems. 
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Our  original  investigation  of  ion-pair  formation  in  Ar*  +  I2  was 
initiated  with  the  expectation  of  testing  the  alkali  analogy  and  probing  for 
possible  effects  caused  by  the  ion-pair  formation  process  occurring  completely 
within  an  ionization  continuum.  In  the  original  publication  (Appendix  A) 
supported  by  this  contract,  we  showed^  the  striking  overall  similarity  between 
the  ion-pair  formation  differential  cross  sections  for  the  reactions 

Ar*  +  I2  ♦  Ar+  +  I2~  (5) 


and 


K  +  I2  *  K+  +  I2"  (6) 

at  center-of-mass  energies  of  ~  25-130  eV.  This  similarity  occurred  despite 
both  the  competing  Penning  channel  {reaction  (1)  of  Appendix  A]  and  the  fact 

it 

that  the  inner-shell  vacancy  in  the  Ar  projectile  implied  multiple  sets  of 
interacting  covalent  and  ion-pair  surfaces  in  contrast  to  the  single  surface 
crossing  in  the  alkali  interaction.  Vfe  explain  this  result  with  reference  to 
the  very  large  Ar-I2  distance  at  which  the  electron  transfer  occurs.  As  noted 
in  Appendix  A,  the  ion-pair  cross  section  is  so  large  that  the  most  important 
contributions  occur  at  impact  parameters  much  larger  than  those  that  would  be 
expected  for  a  Penning  process.  Also,  at  the  large  surface  crossing  distances 
important  to  ion-pair  formation,  the  core  electron  vacancy  should  have  only  a 
very  minor  effect  on  the  potential  shapes  and  surface  coupling  matrix 
elements;  deviations  from  alkali-like  properties  should  only  occur  at  smaller 
distances  where  core  overlap  becomes  significant  (see  Appendix  B). 
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More  quantitative  comparisons,  with  appraisals  of  the  effects  on  the 
distributions  of  reaction  (5)  caused  by  the  competing  Penning  ionization  and 
excitation  transfer  channels,  were  initially  hampered  by  the  lack  of  a  compu¬ 
tational  model.  A  simple  classical  trajectory,  surface-hopping  model  had  been 
used  previously  by  Los  and  coworkers to  treat  ion-pair  formation  in  reaction 
(2).  He  modified  and  extended  their  model  to  include  the  effects  of  depletion 
by  a  competing  continuum  channel.  Our  calculations^  treated  the  continuum 
process  as  an  absorption  width  that  was  assumed  to  be  important  only  along  the 
original  covalent  surface,  not  the  ion-pair  surface.  The  calculations 
demonstrated  the  small  differences  between  reactions  (5)  and  (6),  as  well  as 
the  most  reasonable  explanations  for  the  observed  differences.  This  work  is 
summarized  in  Appendix  B. 

The  overriding  utility  of  the  alkali  analogy  cannot  be  overemphasized, 
because  the  results  for  one  reaction  (Ar*  +  I2)  led  us  to  conclude  that  the 
essential  features  of  the  ion-pair  formation  product  distributions  in  all 
Rg*  +  %2  systems  can  be  predicted  from  comparisons  with  suitably  chosen  M  +  X2 
data.  Since  a  large  body  of  systematic  data  already  exist  for  M  +  X2  inter¬ 
actions,  our  conclusion  vastly  expands  the  impact  of  the  work  on  Ar*  +  I2. 

The  extrapolation  to  other  Rg*  +  X2  systems  was  consistent  with  our  general 
desire  to  achieve  predictive  Insight  from  studying  carefully  chosen  represen¬ 
tative  trlatomic  systems.  It  also  induced  us  to  direct  our  subsequent  effort 
into  collision  systems  with  somewhat  different  properties,  especially  those 
with  lon-pair  formation  cross  sections  that  are  not  so  dominant  relative  to 
other  processes. 

Before  turning  to  other  systems,  we  should  comment  on  the  structure 
observed  (see  Figure  1  of  Appendix  A)  In  the  angular  distribution  for  reaction 
(5).  At  high  energies,  two  peaks  are  evident.  The  "Ionic"  peak  at  large 
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angles  Is  a  rainbow  feature  derived  from  trajectories  that  transfer  to  the 
strongly  attractive  coulombic  (ionic)  surface  as  the  incoming  system  first 
crosses  the  Intersection  region  (see  Figure  2  of  Appendix  A).  The  smaller 
angle  "covalent"  peak  contains  trajectories  that  remain  on  the  covalent  sur¬ 
face  until  they  become  ion-pairs  on  the  second  (outward)  crossing  of  the  two 
surfaces. 

At  lower  energies,  the  dominant  ionic  rainbow  peak  develops  additional 
structure,  a  shoulder  on  its  bright  (low-angle)  side.  In  many  other  recent 
measurements  of  ion-pair  formation  differential  cross  sections,  a  similar 
shoulder  appears  on  the  bright  side  of  the  main  ionic  rainbow  peak.  Although 
this  feature  also  appears  in  trajectory  calculations  on  these  systems,  its 
essential  cause  defied  explanation  until  we  derived  an  analytical  description 
of  the  differential  cross  section  for  a  very  simple  model  two-state  system. 

The  derivation,  an  explanation,  and  a  sample  calculation  for  the  Na  +  I  system 
are  given**  in  Appendix  C. 

Ar*  +  Qj  Ion-Pair  Formation 

The  large  vertical  electron  affinity  (EAy)  of  the  I2  molecule  (and  other 
halogens)  is  the  major  cause  of  the  large  ion-pair  cross  section  because  the 
asymptotic  energy  separation  between  the  covalent  and  coulombic  surfaces  is 

AE  -  IP  -  EAy  (7) 

and  the  surfaces  cross  at  a  distance 

R.  -  1/AE,  (8) 
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where  IP  is  the  ionization  potential  of  the  projectile  and  Eq.  (8)  is  in 
atomic  units;  R£  is  10.6  a.u.  for  Ar*  +  I2<  Significantly  smaller  crossing 
distances  would  imply  smaller  ion-pair  formation  cross  sections.  Oxygen 
molecules  have  a  much  smaller  electron  affinity  than  I2  and,  therefore,  yield 
much  smaller  ion-pair  cross  sections  in  collisions  with  M  or  Rg*.  Hence 
competing  processes  might  be  expected  to  have  more  pronounced  effects  on  the 
observed  distributions. 

We  have,  therefore,  investigated  ion-pair  formation  in  the  Ar*  +  O2 
system. 


Ar*  +  Oj  Ar+  +  02“.  (9) 

Differential  cross  section  measurements  for  this  reaction  are  given  in 
Figure  1.  Our  efforts  on  this  system  were  closely  coordinated  with 
simultaneous  measurements  of  ion-pair  formation  in  the  analogous  reaction 

K  +  02  ♦  K+  +  0^  (10) 

1  O 

by  Kleyn  et  al.  at  FOM.  Additional  data  on  reaction  (10)  were  obtained  by 
13 

Young  et  al.  at  Argonne.  The  available  data  for  reaction  (10)  are  shown  in 
Figure  2. 

The  excitation  energy  of  Ar*  is  slightly  lower  than  the  ionization  energy 
of  02.  Hence  the  continuum  Ar  +  02+  +  e~  is  asymptotically  slightly  above  our 
incoming  channel  (see  Figure  3).  However,  we  still  expected  a  large  cross 
section  for  excitation  transfer  to  the  near-resonant  Rydberg  states  of  02, 

Ar*  +  0j  ♦  Ar  +  02**,  (11) 
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Figure  1.  Experimental  Ar*  +  O2  ion -pair  foraatlon  differential  croea  section 
data  from  our  laboratory.  Each  curve  la  labeled  with  the  energy 
(E0)  of  the  Ar*  beam. 
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Figure  2.  Experimental  K  +  0j  ion-pair  formation  differential  croaa  uetloa 
data  from  PGM  (E.  -  35.0  and  90.5  e?)  and  Ar gonna  (all  otbara). 


a  pseudo-continuum  process,  as  discussed  by  Miller  and  Morgner.*4  In  fact,  we 
have  used  attenuation  techniques  to  show  that  reaction  (9)  is  less  than  >)Z  of 
the  aeasured  total  (20-25  A*)  destruction  cross  section  for  Ar  in  collisions 
with  02. 

Recently,  using  approximate  formulas  of  Miller  and  Morgner,14  we  have 
estimated  the  pseudo-continuum  coupling  due  to  reaction  (11)  and  have  com¬ 
pleted  a  classical  trajectory  study  of  the  pair  of  reactions  (9)  and  (10). 
Although  the  estimated  cross  section  for  reaction  (11)  (~  15-20  k  )  is 
substantially  larger  than  the  cross  sections  (5-7  A  )  for  reactions  (9)  and 
(10),  it  surprisingly  causes  only  a  small  difference  in  the  shape  of  the 
calculated  differential  cross  section.  Since  the  differences  in  shape  between 
our  Ar*  +  02  data  (e.g..  Figure  1,  40  eV)  and  the  K  +  02  data  are  lesa  than 
the  differences  between  the  FCM  data1^  (e.g..  Figure  2,  35  eV)  and  the  Argonne 
data13  (e.g..  Figure  2,  38.2  eV),  one  could  hardly  expect  to  extract  any 
Information  on  the  pseudo-continuum  reaction  from  a  comparison  of  Figures  1 
and  2,  and  the  trajectory  calculations  indeed  confirm  this  observation.  Once 
again,  the  ion-pair  formation  process  in  Ar*  +  Oj  is  not  significantly 
different  from  expectations  based  on  an  analogy  to  alkali  systems.  A 
manuscript  detailing  this  work  will  be  prepared. 

Vibrational  Timing  in  Transient  0," 

The  importance  of  the  ion-pair  surface  is  not,  however,  limited  to  reac¬ 
tions  involving  ion-pair  formation.  In  fact,  the  ion-pair  surface  can  act  as 
a  transient  intermediate  that  profoundly  affects  inelastic  scattering  as  well 
(both  vibrationally  and  electronically).^  Since  the  equilibrium  bond  length 


13 


of  (>2~  is  larger  than  that  of  02,  the  electron  transfer  from  Ar*  to  O2  at  the 
surface  crossing  generates  an  Initially  compressed  (see  Pigure  4)  that 
then  expands  and  begins  to  vibrate  as  the  collision  evolves. 

For  Ar*  +  O2  collisions  in  the  energy  range  that  we  investigate,  the 
vibrational  period  of  the  transient  02~  ion,  formed  after  the  initial 
(incoming)  transfer  to  the  ion-pair  surface,  is  comparable  to  the  collision 
time.  The  path  chosen  by  the  system  at  the  second  (outgoing)  crossing  of  the 
ionic  and  covalent  surfaces  will  depend  profoundly  on  the  state  of  the 
transiently  formed  02”  ion  at  the  Instant  of  that  second  crossing.  This  is 
because  the  vertical  electron  affinity  EAV  changes  dramatically  with  02_  bond 
length  (see  Figure  4),  and  EAy  determines  the  crossing  distance  [from  Eqs.  (7) 
and  (8)]  on  the  outgoing  trajectory,  which  in  turn  profoundly  alters  the 
matrix  element  coupling  the  two  interacting  surfaces. 

Appendix  D  explains  these  effects  in  more  detail,  while  describing15  the 
first  observations  of  large  Internal  excitation  of  neutral  products  mediated 
by  02"  expansion  on  a  transiently  formed  ion-pair  surface.  Similar  results 
have  also  been  reported1 ^  for  alkali -oxygen  collisions. 

Our  efforts  in  the  past  few  years  have  convinced  us  that  these  previously 
neglected  transient  timing  effects,  caused  by  bond  expansions  and  vibrations 
in  fleeting  collision  intermediates,  should  be  Important  to  the  dynamics  of 
many  other  collision  systems.  Appendix  E  reviews^  the  considerable  range  of 
observations  on  transiently  formed  02~  and  mentions  results  for  several  other 
systems.  Of  particular  note  is  the  process  of  covalent  reneutrallzation1^ 

(see  Figures  7  and  8  of  Appendix  E),  where  expansion  of  a  transiently  formed 
12”  induces  a  third  surface  crossing  and  reneutralization  of  Ar+  +  I2“  in  the 
Ar*  +  I2  collision. 
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Figure  4.  Schematic  representation  of  the  Ar*  +  0^  potential  surface  Inter¬ 
secting  the  couloabic  Ar'*’  +•  O^-  surface:  Only  the  lower  surface  is 
shown  [covalent  at  large  Ar-Oj  separation  R,  couloabic  (shadad)  at 
small  R] ,  although  dashed  lines  partially  outline  the  upper 
surface.  A  dashed  trajectory  is  shown  entering  along  the  valley 
from  cop  right  to  bottom  left  (at  fixed  0-0  separation  r),  crossing 
the  surface  Intersection  seam,  Chen  expanding  in  the  r  direction 
and  vibrating  as  it  heads  back  to  large  R  and  a  second  seam 
crossing.  The  upper  insert  is  a  cut  through  the  surface  at  r  »  r# 

( C) ) ,  showing  an  ionic  trajectory  to  ion -pair  products.  The  lower 
Insert  shows,  the  Oj  and  0^  potential  curves;  the  vertical  energy 
separation  between  the  curves  (RAV)  depends  strongly  on  r. 
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Orientational  Rainbows 


For  Ax*  +  I2  collisions  we  have  noted  that  the  peak  at  larger  angles  In 
the  angular  distribution  for  ion-pair  foraation  Is  a  rainbow  associated  with 
"Ionic  trajectories"  that  travel  along  the  couloablc  surface  Inside  Rc.  The 
second  peak  at  lower  angles  is  doainated  by  contributions  froa  the  alternative 
"covalent  trajectories."  Surprisingly,  the  K  +  O2  and  Ar*  +  Oj  trajectory 
calculations  indicate  a  very  large  contribution  to  the  low  angle  peak  froa 
ionic  trajectories.  An  example  is  shown  in  Figure  5,  in  which  the  calculated 
differential  cross  section  is  partitioned  into  its  contributions  froa  the  two 
possible  paths  to  ion-pair  products.  The  new  and  unusual  feature  is  associ¬ 
ated  with  "ionic"  trajectories  having  lapact  paraaeters  slightly  saaller  than 
those  responsible  for  the  principal  rainbow  peak,  where  in  the  analogous  atow- 
atoa  case,  there  is  a  zero  in  the  deflection  function  (with  attractive  and 
repulsive  forces  just  balancing).  For  these  trajectories  in  the  atoa-aolecule 
case,  the  attractive  and  repulsive  forces  are  nearly  balanced,  but  in  general 
the  orientation  of  the  target  aolecule  reaoves  the  balance  and  produces  a 
deflection  function1®  with  a  (non-zero)  minimum  deflection  angle. 

Figure  6  shows  calculated  deflection  functions  for  Ar*  +  02  collisions 
for  two  molecular  orientations.  In  each  case,  the  dashed  portion  and  the 
solid  portion  are  the  covalent  branch  and  the  ionic  branch  of  the  deflection 
function,  respectively.  For  the  light  curves  (a  special  case),  the  deflection 
functions  go  through  zero  at  the  transition  from  attractive  (large  b)  to 
repulsive  (saaller  b)  interactions.  In  the  more  general  orientation  (heavy 
curves),  there  is  a  minimum  in  the  deflection  function;  we  have  named  this 
minimum  the  orientational  rainbow  (OR).  Even  after  averaging  over  all 
orientations,  the  calculations  yield  a  peak  in  the  ion-pair  foraation 


16 


0  100  200 


r  (eV-deg.) 

Figure  5.  Calculated  ion-pair  formation  differential  cross  sections  for  K  + 
Oj  *t  E0  ■  33  eV.  Shown  also  are  the  separated  contribution  for 
the  Ionic  and  covalent  paths,  the  oaln  rainbow  Is  marked  R;  the 
orientational  rainbows  are  denoted  OR.. 
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Figure  6.  Calculated  deflect ion  functions  for  ion -pair  formation  in  Ax*  +  0, 
assuming  two  different  orlentatlona  of  the  Oj  target.  The  covalent 
branch  of  each  curve  is  dashed.  R  indicates  the  classical  rainbow, 
and  CR  locates  the  two  orientational  rainbows  observed  for  one  of 
the  O2  orientations. 


cross  section  consistent  with  this  effect.  In  fact,  a  second  orientational 
rainbow,  associated  with  the  covalent  path.  Is  also  Indicated  by  Figures  5  and 
6. 

Our  trajectory  calculations  for  ion-pair  formation  in  K  +  02  and  •■*  ♦  Oj 
yield  another  surprising  and  Interesting  result:  a  strong  energy  dependence 
In  the  ratio  of  the  two  major  peaks  in  the  differential  cross  section  for  lon- 
pair  production.  As  the  collision  energy  is  lowered,  there  Is  a  dramatic  jump 
In  the  relative  intensity  of  the  low  angle  peak.  This  jump  Is  caused  by  a 
rapid  Increase  In  the  ionic  orientational  rainbow  peak  and  can  be  understood 
by  examining  a  vibrational  timing  diagram  (Figure  7,  cf.  discussion  of  Figure 
4  In  Appendix  D). 

The  impact  parameters  contributing  to  the  orientational  rainbow  are 
confined  to  a  small  range  between  ~  2.3  and  2.6  A.  At  high  energies  (e.g.,  90 
eV),  the  ionic  trajectories  reach  the  second  surface  crossing  after  less  than 
a  single  02~  vibration,  and  the  compressed  02”  of  relatively  small  electron 
affinity  yields  a  significant  reneutralization  and  concomitant  loss  of  the 
Ions  contributing  to  the  OR  feature.  As  the  energy  Is  lowered,  these 
trajectories  evolve  more  slowly  relative  to  the  02”  vibrational  period;  at  low 
enough  energies  (e.g.,  40  eV  for  b  -  2.5  A  in  Figure  7),  the  02“  goes  through 
more  than  one  vibration  before  the  second  surface  crossing.  The  crossing  then 
occurs  with  an  expanded  02~  of  large  EAy,  and  a  much  larger  fraction  of  the 
Initial  ion-pairs  pass  through  the  crossing  without  reneutralizing. 

Figure  8  shows  the  calculated  ion-pair  peak  ratios  for  Ar*  +  02  and 
K  +  02  over  a  beam  energy  range  from  30  eV  to  90  eV.  Also  shown  in  Figure  8 
are  the  comparable  data  for  K  +  O2  (Figure  2)  and  Ar*  +  O2  (Figure  1)  over  a 
slightly  larger  energy  range.  There  Is  evidently  a  very  good  qualitative 
correspondence  between  the  data  and  the  calculation.  The  slight  shift  In  the 
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Figure  7.  Timing  diagram  for  Ax  +  Cw.  The  zero  of  Che  time  scale  represents 
Che  inscant  of  02~  formacion  as  Che  incoming  system  crosses  co  che 
ionic  surface.  Ac  chls  instant,  che  Oj”  begins  vibrating,  causing 
synchronous  oscillations  in  che  vertical  electron  affinity  and  che 
croeaing  distance  Rc';  Re'(t)  is  shown  by  che  heavy  periodic 
curve.  The  light  curves,  marked  first  by  their  energies  (In  eV) 
and  second  by  their  impact  parameters  (in  A),  show  various 
illustrative  trajectories  by  plotting  che  Ar**^'  distance,  R(t). 
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Figure  8.  do/du  sine  peak  ratios:  low  angle  (covalent)  reglon/ionlc  rainbow 
region  plotted  vs.  bean  energy. 

Solid  line:  K  +  O2  calculations 

Dashed  line:  Ar  +  Oj  calculations,  no  contlnuua  depletion 

Triangles:  Ar*  +  Oj  data 

Open  circles:  K  +  Oj  data  (Ar goons) 

Closed  circles:  K  +  Oj  data  (FQK). 
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cricical  transition  energy  indicates  that  the  02~  vibrational  notion  nay  be 
delayed  or  slowed  in  the  field  of  the  positive  ion  relative  to  the  notion 
expected  for  an  isolated  02”. 

The  unusual  effect  shown  by  Figure  8  had  not  been  observed  previously 
and,  although  complicated,  can  be  described  in  a  purely  classical  framework. 
This  process  is  an  energy -dependent,  vlbratlonally  clocked,  orlentatlonally 
derived,  rainbow  feature  in  an  electronically  Inelastic  transition,  and  it 
represents  a  coupling  of  electronic,  vibrational,  rotational,  and  transla¬ 
tional  motions!  Recent  trajectory  calculations  by  G.  Parlant  (unpublished)  on 
more  realistic  potential  energy  surfaces  have  confirmed  our  observations, 
which  are  being  prepared  for  publication. 


He*  +  09 

Appendix  F  describes20  our  development  of  a  modification  of 
neutralizatlon-relonlzatlon  spectroscopy  that  allowed  simultaneous  and 
independent  measurements  of  differential  cross  sections  for  the  ion-pair 
formation  reactions  of  the  two  metastable  states  of  He  with  02, 

He*(21S)  +  02  *  He+  +  02"  (12) 

and 

He*(23S)  +■  02  ♦  He+  +  02_.  (13) 

Charge  transfer  of  He+  in  Rb  y. e.  ted  a  mixture  of  the  two  metastables; 
the  He*(2lS)  formed  directly  in  the  charge  transfer  and  the  He*  (23S)  formed 


only  aftei  'ascade  from  a  He*(2^P)  precursor. ^  After  the  subsequent  lon-pair 
formation,  the  He+  from  reaction  (13)  was  lower  in  energy  than  the  He*  from 
reaction  (12)  by  an  amount  equal  to  that  of  the  He(2^P  ♦  2^S)  transition. 

Since  the  product  energy  distributions  were  reasonably  sharp,  the  two  reac¬ 
tions  were  easily  resolvable.  As  a  checlc,  we  placed  Na  in  the  charge  transfer 
cell  to  produce  and  react  a  beam  of  pure  He*(2^S).  The  results^  were  in 
excellent  agreement  with  the  distributions  for  reaction  (13)  obtained  using  Rb 
charge  transfer  vapor. 

Figure  9  shows  time-of-flight  measurement  of  the  He*  beam  produced  by  Rb 
and  by  Na  charge  transfer.  Figure  10  shows  representative  Re+  product  energy 
distributions  at  various  laboratory  scattering  angles  for  a  mixed  beam  of 
He  (2  S)  and  He  (2  S)  formed  in  Rb.  In  each  spectrum,  the  peak  at  lower 
energy  is  from  reaction  (13). 

The  measured  differential  cross  sections  for  the  two  ion-pair  formation 
reactions  are  shown  over  a  large  energy  range  in  Figures  11  and  12.  An  expla¬ 
nation  for  the  observed  structure  and  the  remarkable  differences  between  the 
two  reactions  is  given  in  Appendix  F.  The  He*(2*S)  +  O2  reaction  shows  strong 
evidence  for  trajectories  Chat  couple  ehree  potential  surfaces  in  a  single 
trajectory  leading  to  ion-pair  products.  These  trajectories  (bdfdc  in  Figure 
3  of  Appendix  F)  were  also  the  first  observed  collisions  in  which  three 
separate  election  jumps  could  be  inferred. 

Attempts  to  model  these  reactions  with  classical  trajectory  calculations 
required  an  extension  of  the  program  to  Include  three  coupled  states.  Figures 
13  and  1A  compare  the  resulting  calculations  with  experimental  differential 
cross  section  data  at  a  beam  energy  of  70  eV.  The  results  are  only  qualita¬ 
tively  comparable.  Since  the  matrix  elements  coupling  the  surfaces  are 
comparable  in  size  to  the  splitting  between  the  two  Incoming  covalent  surfaces 
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Figure  9.  Tiea-of-flight  aeaaureaants  of  Be*  beaa  coapoeltlon  from  He4*  charge 
transfer  in  Rb  and  Na.  The  two  states  in  Rb  arc  not  resolved  here, 
but  have  been  in  M.  Barat's  laboratory. 
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Figaro  10.  Ioa-palr  formation  doablo  differential  croca  sactlooo  for  Bo*  +  Oj 
colli alone  at  a  boos  anargy  of  70.6  oV  (2*3  scoeo)  and  70.2  «V  (23S 
•cata).  AC  aaeh  laboratory  scat taring  angla  (aarfcad  at  right  of 
aaeh  apactroa),  tba  lower  anargy  paafc  la  from  raactloa  (13),  tha 
hlgfaor  paafc  froa  raaetlon  (12).  Spaetra  art  all  noraallzed  to  tha 
aaa a  lntanalty  at  thalr  paafca;  relative  intaaaleias  can  bo  obtalnad 
froa  Flgusoo  13  and  14. 
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Figure  11.  Differential  cross  section  results  for  reaction  (13)  plotted  vs. 

the  reduced  scattering  angle  t  (  -EQ9).  Beau  energies  (E0)  are 
Indicated  for  each  curve* 
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Figura  12.  Data  for  raacclon  (12).  Saa  Lagand  for  Flgura  11. 
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figure  13.  A  comparison  of  calculation  and  experiment  for  reaction  (13)  at 
-  70.2  eV. 


Figure  14.  A  comparison  of  calculation  and  experiment  for  reaction  (12)  at 


(shown  In  Figure  3  of  Appendix  F),  our  use  of  a  calculation  that  assumes 
isolated  crossings  and  calculates  deflections  using  diabatlc  potentials  is 
obviously  Inadequate.  In  particular,  if  we  modified  our  trajectory  calcula¬ 
tions  to  use  adiabatic  potential  surfaces,  we  anticipate  that  the  central  peak 
in  Figure  14  would  shift  to  significantly  lower  angles.  If  we  also  increased 
the  repulsion  in  the  various  calculated  potential  surfaces,  we  could  then 
probably  achieve  reasonable  agreement  in  angular  position  between  the 
calculations  and  the  experiments. 

However,  the  disagreement  in  relative  peak  Intensities  is  a  such  more 
difficult  problem  because  the  Individual  Intensities  can  depend  strongly  on 
their  respective  depletion  by  competing  continuum  channels. 

He*  +  O2  ♦  He  +■  (>2+  +  e~.  (14) 

*  1 
In  this  case,  the  large  He  excitation  energies  (19.8  and  20.6  eV  for  the  2  S 

and  2*S  states,  respectively)  cause  the  interacting  discrete  potential  sur¬ 
faces  (shown  schematically  in  Figure  3  of  Appendix  F)  to  be  nearly  resonant 
with  several  different  02+  electronic  states  in  the  continuum  channel 
[Equation  (14)].  This  would  suggest  a  large,  perhaps  dominant,  role  of  the 
Penning  process  in  these  systems.  In  fact,  comparison  of  our  attenuation 
measurement  of  total  He*  destruction  with  our  measured  yields  of  fast  positive 
(He+)  and  slow  positive  ((>2+,  0+)  and  negative  (O2”,  0~,  e“)  products 
indicates  that  the  ion-pair  channel  is  only  ~  102  of  the  total  destruction 
cross  section  at  collision  energies  near  100  eV.  The  Penning  channel  domi¬ 
nates  in  this  energy  range.  Results  at  higher  energies  from  another  labora¬ 
tory  are  in  agreement  with  this  conclusion. 22 


In  addition,  for  a  recent  thermal  energy  investigation23  of  reaction 
(14),  product  electron  energy  distributions  were  interpreted  by  assuming  a 
significant  continuum  depletion  of  reactants  while  they  were  on  the  ionic  He* 
+  02*  surface.  If  this  is  true,  our  modeling  in  the  trajectory  calculations 
of  the  depletion  by  continuum  processes  would  need  separate  estimates  of  the 
continuum  widths  for  each  of  the  three  surfaces  explored  in  the  ion-pair 
process.  Although  we  could  attempt  this  by,  for  example,  treating  the  contin¬ 
uum  equivalently  for  all  three  surfaces,  we  find  this  hard  to  justify  in  light 
of  our  earlier  assumption^  in  the  Ar*  +  I2  system  that  there  was  no  continuum 

coupling  from  the  ion-pair  surface.  Clearly,  additional  experimental  informa¬ 
tion  on  the  continuum  interaction  is  sorely  needed  for  any  attempt  to  model 
the  rich  and  complex  He*  +  Oj  interaction. 

Ion-pair  formation  has  also  been  investigated  briefly  for  the  interac¬ 
tions  of  He*(23S)  and  He*(21S)  with  NO,  NO2,  and  N2.  The  neutrallzatlon- 
relonlzatlon  procedure  appears  from  these  reactions  to  be  of  general  utility 
in  studying  He  (2  S)  and  He  (2*S)  simultaneously,  but  separately. 

The  qualitative  results  for  NO  targets  are  similar  to  those  for  O2 ,  as 
was  anticipated  from  the  predicted  similarities  of  the  shapes  of  the 
interacting  potential  surfaces.  With  N2  target  gas,  the  cross  section  is 
quite  small  because  the  negative  electron  affinity  of  N2  implies  a  very  small 
crossing  distance  Rc  between  the  incoming  covalent  surface  and  the  product 
ion-pair  surface  He+  +  N2~.  In  contrast,  NO2  has  a  large  positive  electron 
affinity  and  a  very  large  cross  section  for  ion-pair  formation. 

In  analogy  with  O2  targets,  we  expect  that  the  Penning  processes  are  more 
important  than  ion-pair  formation  in  all  of  these  systems.  For  this  reason, 


detailed  Investigation  of  ion-pair  formation  in  these  systems  was  abandoned  in 
favor  of  an  attempt  to  develop  techniques  for  examining  the  continuum 
processes  directly. 


Direct  Continuum  Channel  Measurements 

In  this  work,  we  have  made  several  significant  contributions  to  the 
understanding  of  the  ion-pair  formation  process  and  the  Inelastic  processes 
that  involve  the  ion-pair  surface  as  a  transient  intermediate.  In  contrast, 
the  competing  continuum  channel  (Penning  ionization)  has  only  grudgingly 
yielded  a  small  amount  of  information. 

The  differential  cross  section  comparisons  between  Ar*  and  K  for  I2  and 
O2  targets  showed  only  small  differences,  and  the  calculations  show  that  these 
distributions  are  not  very  sensitive  to  variations  in  the  competing  channels. 
For  the  He*  interactions,  we  have  demonstrated  the  ability  to  control  the  beam 
composition  and  investigate  He*(2*S)  and  He*(2^S)  interactions  Independently. 
Although  this  gives  us  the  flexibility  for  each  molecular  target  to  investi¬ 
gate  reactions  on  a  coupled  set  of  interacting  potential  surfaces  starting 
from  two  separate  and  well-characterized  input  channels,  it  still  does  not 
yield  any  information  on  the  Penning  channel  that  we  infer  to  be  the  dominant 
inelastic  process. 

Any  attempt  to  measure  differential  cross  sections  for  either  the  atomic 
or  ionic  products  of  the  Penning  reaction  is  hampered  by  a  major  experimental 
difficulty.  Consider  reaction  (14)  as  an  example.  If  we  try  to  detect  the 
scattered  He  atom  product  (possible  at  energies  above  ~  100  eV  where  the 
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secondary  electron  ejection  coefficient  be cooes  reasonably  large),  we  mist 
contend  with  a  very  large  background  of  elastically  and  inelastlcally 
scattered  He*, 

He*  +  O2  He*  +  02^*\  (15) 

Interfering  with,  and  probably  swamping,  the  Intended  measurement.  In 
contrast,  to  determine  differential  cross  sections  for  the  ionic  product(s) 
that  are  produced  with  near-thermal  energies,  we  must  expend  tremendous 
efforts  In  eliminating  contact  potentials,  keeping  surfaces  clean  and 
uncharged,  and  transporting  and  analyzing  the  products  without  distorting  the 
nascent  distributions.  Only  one  group,  working  on  thermal  energy  Penning 
reactions,  has  achieved  any  success^  measuring  the  distributions  of  these 
slow  Ions.  We  have  developed  an  alternative  approach,  suitable  at  higher 
energies,  that  separates  the  He  from  the  He*. 

The  key  to  our  technique  is  the  observation  that  the  fundamental  marker 
for  the  continuum  process  Is  the  formation  of  a  positive  ion  from  the  target 
species.  Hence,  If  this  ion  can  be  extracted  from  the  collision  zone  and 
counted.  Its  appearance  can  be  used  to  validate  the  detection  of  He  atoms  from 
reaction  (14)  and  discriminate  against  He*  from  reaction  (15).  The  two 
products  of  reaction  (14)  are  detected  in  delayed  coincidence  on  two  separate 
detectors. 

A  fast  Be*(23S)  beam,  produced  by  charge  transfer  of  He+  in  Na  vapor, 
crosses  a  beam  of  02  effusing  from  a  multi capillary  array.  Grids  surrounding 
the  interaction  zone  are  used  to  extract  slow  positive  product  ions  to  a 
Channeltron  detector  near  the  intersection  zone.  A  second  Channeltron  can  be 
positioned  at  different  angles  around  the  intersection  zone  to  detect 
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scattered  fast  neutral  products  (He  and  He*).  The  signal  from  the  first 
Channeltron,  after  amplification,  starts  a  time-to-amplltude  converter 
(TAC).  The  stop  Input  of  the  TAC  receives  amplified  pulses  from  the  second 
Channeltron.  If  a  stop  pulse  Is  received  within  the  time  window  of  the  TAC, 
an  output  pulse  (amplitude  proportional  to  the  measured  time  difference)  from 
the  TAC  passes  through  an  analog-to-dlgltal  converter  to  a  multichannel 
analyzer  (MCA)  used  as  a  pulse-height  analyzer.  In  this  way,  the  MCA  produces 
a  spectrum  of  the  time-correlation  of  delayed  coincidences  between  the  two 
detectors. 

(Xir  grid  geometry  and  acceleration  techniques  for  extracting  the  slow 
positive  ions  are  designed  to  minimize  the  time  spread  (for  a  given  ion  mass) 
associated  with  the  product  ion  and  also  to  compensate  for  the  time  spreads 
due  to  the  variation  in  location  of  product  formation  across  the  beam  inter¬ 
section  zone.  Then  the  measured  time  differences  for  reaction  (14)  can  be 
directly  correlated  with  the  flight  times,  and  hence  the  velocities,  of  the 
coincident  He  atoms.  He*  from  process  (15)  is  not  correlated  in  time  with 
product  ions  and  therefore  contributes  only  a  uniform  background  rate  to  the 
coincidence  spectra.  The  very  large  signal  (especially  at  small  scattering 
angles)  from  elastic  scattering  is  discriminated  against  by  spreading  it 
uniformly  over  the  time  axis;  most  of  the  He*  does  not  even  appear  within  the 
small  TAC  time  windows  selected  by  the  initiating  ion-detection  pulse. 

Figure  15  is  a  typical  coincidence  spectrum,  taken  directly  from  the  MCA 

it  o 

without  analysis  or  smoothing,  for  the  He  (25)  +  (>2  Penning  ionization 
reaction.  The  beam  energy  was  296  eV,  and  the  He  product  was  detected  at  a 
laboratory  scattering  angle  of  1.3°.  Rising  above  a  uniform  background  from 
uncorrelated  processes  are  two  peaks  identified  with  the  Penning  processes, 
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Figure  15.  Experimental  delayed  coincidence  spectrua  between  extracted  slow 
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positive  Ions  and  scattered  fast  neutrals  In  Be  (2  S)  +  Oj 
collisions. 
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He*(23S)  +  02  -*■  He  +  02+  +  e~,  (16) 

and 

He*(23S)  +  02  >  He  +  0+  +  0  +  e“.  (17) 

The  lighter  0+  ions,  from  Penning  formation  of  dissociative  or  predissociated 
02+  states,  are  extracted  more  rapidly  than  the  02+  ions  and  therefore  imply  a 
larger  time  delay  before  detection  of  correlated  He  atoms.  It  is  evident  from 
Figure  15  that  we  are  able  to  adjust  the  extraction  and  timing  conditions  so 
that  we  can  examine  reactions  (16)  and  (17)  simultaneously. 

Coincidence  spectra  for  various  scattering  angles  of  the  He  product  are 
given  in  Figure  16.  The  simple,  single-peak,  timing  spectra  evident  at  small 
angles  for  reactions  (16)  and  (17)  become  more  complex  as  the  scattering  angle 
increases.  For  both  reactions,  approximate  exothermicity  (0)  values  are  indi¬ 
cated  above  the  spectra.  Although  the  main  peak  ia  nearly  thermoneutral  (with 
the  electron  carrying  away  the  excess  energy  separating  the  initial  and  final 
potential  surfaces),  the  second  feature  indicates  an  18-20  eV  conversion  of 
translational  energy  either  to  Internal  energy  or  to  an  Increase  in  the 

product  electron  energy.  Very  little  information  now  exists  concerning  02+ 

£  1 

states  18-20  eV  above  the  initial  He  (2  S)  level,  and  no  obvious  explanation 
can  be  found  for  the  Importance  of  this  peak. 

Figure  17  gives  reduced  p  [  ■  0sin9  ~  J  vs.  t  [  »  8Egj  plots  of  the 
angular  distributions  of  the  He  atoms  from  reactions  (16)  and  (17).  Each 
point  is  an  integral  over  the  product  velocity  distribution  in  the  tlme-of- 
flight  spectrum  at  the  appropriate  angle.  Clearly,  the  dissociative  channel, 
reaction  (17)  is  a  significant  fraction  of  the  total  Penning  channel;  this 
result  is  not  surprising  when  one  observes  that  the  He  (2  S)  Internal  energy 
is  greater  than  the  energy  required  to  dlssoclatlvely  ionize  02. 
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Plgure  16.  Experimental  delayed  coincidence  for  Be  (2JS)  +  O2  colllalons.  The 
vertical  line  dividing  the  spectra  separatee  the  contributions  from 
reactions  (16)  and  (17).  Exothermlclty  values  (Q)  are  Indicated  at 
the  cop  for  each  reaction.  The  laboratory  scattering  angle  for 
each  spectrum  Is  shown  at  the  right.  Approximate  relative 
Intensities  are  shown  by  the  scaling  factors  for  p  on  the  left. 
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Figure  17.  Reduced  (p— r)  differencial  croea  section  data  for  reactions  (16) 
and  (17)  at  SQ  ■  296  eV. 


The  low  angle  thresholds  for  the  Penning  processes  are  at  smaller  angles 
than  those  investigated.  Measurements  at  these  smaller  angles  were  Impossible 
because  the  unscattered  Re*  beam  overwhelmed  the  He  detector.  A  better 
colllmailon  of  the  He*  beam  would  allow  measurements  to  be  made  closer  to  the 
beam.  The  decrease  in  beam  Intensity  accompanying  a  tighter  colliraation  would 
present  no  real  problem  because  the  coincidence  signal  levels  are  very  large 
at  small  angles.  Par  example.  Figure  15  was  generated  in  ~  6  minutes;  satis¬ 
factory  spectra  at  half  that  angle  can  now  be  measured  in  ~  1  minute.  In 
fact,  since  we  attempt  to  collect  a  large  fraction  of  the  product  ions,  we 
have  already  been  required  in  these  measurements  to  reduce  the  beam  intensity 
substantially  to  avoid  very  high  count  rates  on  the  ion  detector.  The 
efficient  collection  of  the  product  ions  makes  data  accumulation  much  faster 
in  these  experiments  than  in  most  other  two-particle  coincidence  measurement 
techniques. 

Our  results  are  the  first  differential  cross  section  measurements  for 
Penning  ionization  at  energies  above  the  thermal  range.  They  are  to  our  know¬ 
ledge  the  only  double-differential  cross  section  measurements  at  any  energy 
for  an  ionization  continuum  channel  in  a  neutral  trlatomic  system.  Although 
these  meausrements  are  preliminary  and  should  be  extended  to  smaller  and 
larger  scattering  angles,  they  indicate  the  success  of  our  newly  developed 
technique  and  demonstrate  its  unique  capabilities  for  rapid  and  direct 
examination  of  the  Penning  channel  in  collisions  at  moderate  energies. 
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Differential  ion-pair  production  cron  sections  for  Ar*-lj  collisions  at  center  of  mass  energies  of  25-133  eV  are  quite 
similar  to  analogous  alkali-  Ij  results  despite  the  addition  of  competing  channels  (Penning  ionization  and  excitation  trans¬ 
fer)  and  possible  multiple  electronic  surface  crossings. 


1.  Introduction 

The  formation  of  ion  pairs  from  the  collision  of 
fast  atoms  with  various  target  molecules  has  been  in¬ 
vestigated  extensively  over  the  past  several  years  [1,2]. 
Detailed  differential  cross  section  information  only  ex¬ 
ists,  however,  for  a  few  prototype  systems  [3- 10], all 
involving  collision  partners  in  their  lowest  electronic 
states.  Here  we  describe  the  first  differential  cross  sec¬ 
tion  measurements  of  collisional  ion-pair  formation 
for  electronically  excited  reactants,  a  study  of  the  re¬ 
action  of  metastable  argon  (Ar*)  with  I2.  This  system 
is  unique  in  that  both  the  incoming  covalent  potential 
surfaces  Ar*  +  I2  and  the  product  ionic  surfaces  Ar+  + 
I2  are  imbedded  in  the  continuum  Ar  +  I2  +  e” . 
Hence,  Penning  ionization, 

Ar*  +  I2  Ar  +  I2  +  e~  ,  (1) 

can  compete  with  the  ion-pair  channel 

Ar*  +  I2  -  Ar*  +  Ij  .  (2) 


2.  Experimental 

The  apparatus  is  a  modification  of  one  previously 
described  [II]  and  will  be  discussed  only  briefly.  An 
Ar*  parent  ion  beam  produced  in  a  discharge  source 

*  Supported  by  the  Office  of  Naval  Research. 

*  Visiting  Scientist  1977-78,  on  leave  from  University  of 
Western  Ontario. 

192 


is  converted  to  Ar*  by  near-resonant  charge  transfer  in 
Rb  vapor  [12].  The  product  fast  neutral  beam,  after 
decay  of  radiative  excited  states,  is  a  mixture  of  the 
two  metastable  states  3P0  and  3P2  and  the  [S0  ground 
state.  After  deflection  of  the  unneutralized  ions  this 
neutral  beam  enters  a  collision  cell  of  length  1.3  cm 
filled  with  I,  vapor  at  300  K  and  at  low  density 
(<1013/cm*).  Product  ions  are  energy  analyzed  and 
counted  as  a  function  of  scattering  angle  by  a  1 27s 
analyzer  and  channeltron  detector  that  rotate  around 
the  collision  cell. 

The  energy  analyzer  can  be  floated  at  various  po¬ 
tentials  to  analyze  parent  or  product  ions  with  high  or 
low  resolution  and  to  determine  the  separate  contri¬ 
butions  to  any  measured  energy  width  due  to  the  ion 
energy  distribution  and  the  energy  analyzer  transmis¬ 
sion  function.  The  en-rgy  analyzer  has  a  measured 
resolution  (fwhm)  of  2.8%  of  the  analysis  energy.  A 
typical  value  for  the  energy  spread  (fwhm)  of  the  par¬ 
ent  ion  beam  is  0.7  eV. 

The  energy  of  the  Ar*  beam  differs  from  that  of 
the  parent  Ar*  by  a  contact  potential  due  to  the  Rb- 
coated  charge  transfer  oven  and  a  small  energy  defect 
in  the  charge  transfer.  We  determine  the  Ar*  beam  en¬ 
ergy  by  a  pulsed  time-of-flight  technique  [11]  using 
another  channeltron  detector.  The  measured  angular 
width  of  the  cylindrically  collimated  Ar*  beam  is  typ¬ 
ically  0.5°  fwhm  in  the  plane  of  motion  of  the  detec¬ 
tor  and  is  probably  of  similar  magnitude  perpendicular 
to  this  plane.  The  estimated  laboratory  angular  accep¬ 
tance  of  the  ion  detector  is  ±0. 1°  full  width  (fw)  in 
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the  plane  of  its  motion  and  ±0.4°  (fw)  perpendicular 
to  that  plane. 

3.  Results 

The  differential  cross  section  for  ion-pair  formation 
is  ptesented  in  fig.  1  for  four  center-of-mass  (c.m.)  col¬ 
lision  energies  Ec m  .  For  each  curve  the  laboratory  an¬ 
gular  distributions  weighted  by  sin  8  are  plotted  ver¬ 
sus  the  reduced  scattering  angle  r  “  E^B.  These  an¬ 
gular  distributions  were  measured  with  the  energy  anal¬ 
yzer  adjusted  to  transmit  the  peak  of  the  product  en¬ 
ergy  distribution  at  each  scattering  angle.  The  observ¬ 
ed  scattering  was  confined  to  low  laboratory  angles 
(<  10°),  and  the  measured  product  energy  profiles  ap¬ 
peared  to  be  almost  independent  of  scattering  angle. 


Fig.  1.  Angular  distributions  of  product  Ar*  lorn  from  colli¬ 
sion  of  Ar*  with  lj  at  four  cantsr-of-nuss  tnargias  £c,m.. 
The  measured  differential  cross  sections  have  been  weighted 
by  sin  and  plotted  afainst  the  reduced  angle  r  ■ 
£c.m.*c.m.-  Successive  curves  have  been  displaced  vertically 
by  0.5  units  for  clarity. 


At  several  widely  spaced  angles,  the  entire  product  en¬ 
ergy  distribution  was  measured.  In  every  case,  a  single 
narrow  peak  was  observed  with  a  width  (fwhm)  of 
£0.8  eV.  The  measured  endothermicities,  A£,  were 
2.6  ±  0.4  eV.  Since  the  product  energy  distributions 
did  not  vary  significantly  with  scattering  angle,  the 
measured  angular  distributions  at  fixed  A E  plotted  in 
fig.  1  are  an  accurate  representation  of  the  total  prod¬ 
uct  ion  angular  distributions.  However,  predominant 
small  angle  scattering  of  the  Ar*  product  ions  com¬ 
bined  with  the  apparatus  angular  resolution  of  >0.5° 
fwhm  implies  a  significant  amount  of  angular  smear¬ 
ing  at  angles  below  <*>  1°.  This  smearing  is  quite  signi¬ 
ficant  at  the  highest  energy  studied. 


4.  Discussion 

A  useful  framework  for  discussing  reactive  en¬ 
counters  involving  rare  gas  metastable  (Rg*)  projec¬ 
tiles  is  based  on  the  observation  [13—19]  of  the  sim¬ 
ilarities  between  Rg*  atoms  and  alkali  atoms  M  due 
to  the  importance  in  both  species  of  the  single  loose¬ 
ly  bound  (4— S  eV)  outer  s  electron.  The  availability 
of  a  large  body  of  detailed  collision  data,  both  exper¬ 
imental  and  theoretical,  involving  alkali  atoms  [1, 
20-22]  often  provides  a  starting  point  for  predicting 
important  features  of  Rg*  interactions. 

For  the  venerable  alkali  reaction  with  a  halogen 
molecule  (X2) 

M  ♦  X2  -  MX  +  X  ,  (3) 

it  it  well  established  [20]  that  the  reaction  dynamics 
at  thermal  energies  are  dominated  by  a  single  poten¬ 
tial  surface  crossing,  at  a  large  intemuclear  separation, 
between  relatively  flat  M  +  X2  diabatic  surface  and 
an  attractive  coulombic  M*  +  X2  surface.  At  the  cross¬ 
ing  distance,  rc,  an  electron  is  said  to  jump  from  M  to 
X2,  forming  X2  on  its  repulsive  inner  wall  in  a  com¬ 
pressed  configuration  that  dissociates  in  the  field  of 
the  positive  ion  as  the  ionic  MX  bond  forms. 

At  collision  energies  significantly  above  thermal, 
dissociation  along  the  coulombic  surface  becomes  en¬ 
ergetically  accessible  and  can  yield  either  M+  +  X~  * 

X  or  M+  ♦  X2  products.  This  ion-pair  formation  chan¬ 
nel  [1,2]  can  still  be  most  simply  understood  in  terms 
of  the  effect  of  crossinp  between  the  two  important 
potential  surfaces,  if  proper  account  is  taken  of  X2 
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bond-stretching  during  the  collision  [23,10]. 

The  corresponding  ion-pair  formation  reactions  of 
Rg*  with  X2  are  similarly  influenced  by  intersections 
of  covalent  Rg*  +  X2  and  coulombic  Rg*  +  X2  sur¬ 
faces  whose  shapes  and  properties  at  large  r  are  very 
similar  to  the  relevant  M-X2  surfaces.  However,  the 
inner-shell  vacancy  in  Rg*  gives  it  unique  properties 
and  interaction  mechanisms  potentially  more  com¬ 
plex  than  those  of  alkali  atoms.  First  there  are  com¬ 
peting  channels  such  as  Penning  ionization,  reaction 
(1),  and  electronic  energy  transfer  (e.g.,  Ar*  + 12  •* 

Ar  + 12)  which  may  or  may  not  perturb  the  ion-pair 
results,  depending  on  the  relative  importance  of  the 
mechanisms.  Second,  the  Ar*  beam  contains  an  un¬ 
known  mixture  of  two  metastable  states  (3P2  and 
3P0)  separated  by  0.17  eV  [24]  with  symmetry  and 
degeneracy  properties  quite  distinct  from  the  M  atom 
(2S).  There  is  a  similar  distinction  between  the  Ar* 
ion  with  two  states  split  by  0.17  eV  (2P3^2  and  2P|/2) 
and  the  M+  ion  ( 1 S).  Although  the  expected  low  prob¬ 
ability  of  core  rearrangement  interactions  strongly 
suggests  that  electron  transfer  will  involve  predomi¬ 
nantly  the  Ar*(3P2)  -*  Ar*(2P3/2)  and  Ar*(3P0)  -*• 
Ar+(2Pi/i)  transitions,  an  additional  complication 
arises  from  the  various  distinct  potential  surfaces  that 
derive  from  interactions  of  each  of  these  four  atomic 
states  with  the  diatomic  (I2  or  I2  ).  All  of  these  cova¬ 
lent  and  coulombic  surfaces  may  not  deviate  signifi¬ 
cantly  in  shape  from  the  corresponding  alkali  sur¬ 
faces,  except  at  small  r,  but  the  transition  rates  at  the 
surface  crossings  (at  large  r)  could  be  quite  different 
due  to  either  the  altered  symmetry  properties  or  the 
interactions  between  a  number  of  closely  spaced  sur¬ 
faces.  Thus,  it  would  not  be  surprising  to  observe 
some  significant  deviations  from  analogous  alkali  re¬ 
sults  in  the  features  of  reaction  (2). 

If  a  value  of  2.6  eV  is  used  (25]  for  the  adiabatic 
electron  affinity  of  I2  a  minimum  endothermicity 
A E  of  * 1.6  eV  is  necessary  to  produce  ion  pairs  with 
the  product  IJ  in  the  lowest  vibrational  level  of  its 
ground  electronic  state  (22J).  Higher  A E  values 
would  correspond  to  production  of  l2  with  internal 
excitation.  A  translational  energy  loss  of  *2.7  eV  is 
needed  to  produce  Ar+  + 1“  (*S)  +  l(2P3^ 2)  products, 
and  A E  *  3.6  eV  is  required  [24]  before  the  Ar*  + 

I-  + 1*  (2Pl/2)  final  state  can  be  reached.  If  the  I2 
internal  energy  is  equal  to  the  energy  generated  by  a 
vertical  transition  from  12  to  I2  at  the  12  equilibrium 
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bond  length  (2.67  A),  the  predicted  A E  value  would 
be  2.S  eV  using  the  appropriate  (vertical)  electron  af¬ 
finity  EAy  of  1.7  eV  determined  by  Hubers  et  al.  [26]. 

The  measured  energy  loss  of  2.6  ±  0.4  eV  implies 
a  near  vertical  transition  (12  -►  12  )  at  the  crossing  of 
the  Ar*  + 12  covalent  surface  with  Ar*  + 12  (2SJ).  It 
also  implies  that  there  is  little  likelihood  of  any  ex¬ 
cited  I*(2Pj/2)  product.  The  majority  of  collisions 
must  therefore  yield  very  highly  vibrationally  ex¬ 
cited  I2  in  the  2S*  ground  electronic  state,  but  the 
range  of  measured  exothermicities  allows  for  the 
production  of  some  dissociated  I~  +  I  products  as 
well  [26]. 

In  the  simplest  approximation  the  crossing  be¬ 
tween  the  “coulombic”  ionic  and  “flat”  covalent  sur¬ 
faces  occurs  at  an  Ar*-I2  distance 

rc  *  14.4  A/(IP  -  EA) , 

where  IP  *  4.21  eV  (the  iohization  potential  of  Ar*), 
and  EA  is  the  electron  affinity  of  I2.  Using  EA,  * 

1.7  eV,  rc  is  approximately  S.7  A;  at  such  large  r  val¬ 
ues,  a  coulombic  approximation  for  the  AR*  +  12  in¬ 
teraction  and  an  assumption  of  negligible  interaction 
between  Ai*  +■  12  ate  both  reasonable.  As  indicated 
schematically  in  fig.  2,  there  are  two  distinct  paths 
(a  and  b)  leading  to  ion  products.  Since  path  b  (the 
ionic  path)  involves  significantly  more  attractive  in¬ 
teraction  than  path  a  (the  covalent  path),  it  should 
produce  Ar*  scattering  at  larger  angles.  At  the  higher 
energies  shown  in  fig.  1,  there  are  indeed  two  major 
peaks  in  the  angular  distribution.  The  r  values  of 
these  peaks  are  very  similar  to  those  found  [7,8]  in 
K  + 12,  the  most  closely  related  M  +  I2  system  inves¬ 
tigated.  As  the  collision  energy  is  lowered,  the  ionic 
peak  in  fig.  1  broadens,  develops  secondary  structure, 
and  increases  in  intensity  relative  to  the  accompanying 
covalent  peak.  These  results  are  again  strikingly  simi¬ 
lar  to  the  published  results  for  K  + 12,  and  the  inter¬ 
pretation  given  below  is  quite  analogous  to  recent  con¬ 
clusions  of  Los  and  co-workers  [2,7,23]  for  M  +  X2. 

A  simple  Landau-Zener  description  of  the  crossing 
region  incorrectly  predicts  that  the  two  major  peaks 
have  comparable  total  intensities.  The  reason  that  the 
ionic  path  has  a  higher  yield  of  product  ions  can  be 
understood  by  examining  the  nature  of  trajectories 
that  explore  the  ionic  surface  inside  the  initial  crossing 
r  =  5.7  A.  As  indicated  in  the  simplified  picture  in  fig. 
2,  the  electron  jump  occurs  to  an  I2  molecule  near 
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Fig.  2.  Sections  through  an  idealized  Ar-I-I  potential  sur¬ 
face.  Upper:  potential  as  a  function  of  Ar-Ij  distance  at  an 
I- 1  distance  (2.67A)  corresponding  to  the  equilibrium  Ij 
bond  separation.  This  is  a  one-dimensional  representation  of 
the  ion-pair  formation  channel  and  shows  two  different 
paths  through  the  crossing  at  S.7  A  leading  to  ion-pair  prod¬ 
ucts.  Lower:  potential  as  a  function  of  I-I  distance  at  Ar-Ij 
distances  of  ••  and  S.7  A.  The  coulombic  Ar  *  *  If  curve 
drops  in  energy  as  'Ar-li  decreases;  at  r/tt-Ij  ■  5-7  A, 
transfer  from  I2  to  If  can  occur,  producing  the  compressed  If. 

its  equilibrium  intemuclear  distance,  forming  an  I2 
molecule  on  its  repulsive  inner  wall.  As  the  Ar+  con¬ 
tinues  to  approach,  reaches  the  turning  point,  and 
then  recedes  from  the  Ij ,  the  12  bond  is  stretching. 

At  the  second  crossing  of  the  ionic  and  covalent  sur¬ 
faces,  the  I2  bond  is  larger,  implying  a  larger  effec¬ 
tive  electron  affinity,  a  larger  crossing  diatance,  and  a 
more  diabatic  crossing  [27],  At  high  collision  veloci¬ 
ties,  the  bond  length  changes  little  between  the  two 
crossings;  as  the  collision  velocity  is  lowered,  the  first 
crossing  becomes  more  and  more  adiabatic  and  bond 
stretching  yields  larger  diabatic  probabilities  (and 
hence  enhanced  ion-pair  formation)  on  the  return 
passage. 

The  additional  structure  observed  at  low  energies 
in  the  ionic  channel  is  more  difficult  to  interpret  [28]. 


A  recent  model  calculation  of  Aten  and  Los  [10]  for 
K  +  Br2  suggests  that  this  structure  is  at  least  partial¬ 
ly  due  to  prestretching  of  the  X2  bond  before  the  re¬ 
actants  reach  the  crossing  seam  at  rc. 

There  is  no  clear  evidence  of  any  significant  observ¬ 
able  effect  on  the  scattering  patterns  in  the  ion-pair 
chancel  (2)  by  the  competing  Penning  or  excitation 
transfer  channels.  This  is  evident  from  the  similarity  of 
the  Ar*+ 12  scattering  to  both  theK  +  I2  results  and 
the  most  recent  model  calculations  [10,28]  forM-F^. 
Further  evidence  is  provided  by  our  preliminary  ob¬ 
servation  that  the  total  destruction  cross  sections  of 
At*  in  I2  are  near  100  A2,  well  above  the  expected 
Penning  cross  sections  for  this  energy  range  [29,30]. 
Basically  the  most  important  contributions  to  the 
charge  transfer  reaction  occur  at  such  large  impact 
parameters  that  they  are  not  significantly  altered  by 
processes  occurring  only  at  shorter  range. 
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Comparison  of  ion  pair  formation  in  the  systems  Ar*+I2 
and  K+l2 

A.  P.  Hickman  and  Keith  T.  Giflen 

Molecular  Physics  Laboratory.  SKI  International  Menlo  Park.  California  94025 
(Received  1  April  1980;  accepted  30  June  1980) 

A  simple  model  that  has  been  used  extensively  by  Los  and  co-workers  to  treat  ion  pair  formation  in  collisions 
of  alkali  atoms  with  diatomic  molecules  is  extended  to  include  continuum  coupling  via  s  competing  Penning 
ionization  channel.  This  extended  model  is  then  used  to  calculate  the  differential  cross  sections  for  ion  pair 
formation  for  the  system  Ar'+I,  over  the  energy  range  28-154  eV  and  to  compare  with  a  previous  treatment 
of  K  +  Ij.  In  the  absence  of  significant  competition  from  continuum  processes,  Ar‘  is  expected  to  behave 
in  a  manner  similar  to  K,  since  the  active  electron  is  an  unpaired  4s  electron  in  both  cases.  We  perform 
model  calculations  for  Ar‘+I,  to  investigate  the  effects  of  varying  the  potential  curves  and  charge 
exchange  matrix  elements  and  of  including  a  continuum  coupling  function  r(R)  Comparison  with 
previous  calculations  for  K  +  I,  suggests  increased  repulsion  on  the  Ar‘-f,  surfaces  relative  to  those  of 
K-I,.  The  competing  mechanisms  of  excitation  transfer  and  Penning  ionization  may  have  a  small  effect 
upon  the  ion  pair  angular  distributions. 


I.  INTRODUCTION 

The  study  of  a  pair  of  generally  similar  reactions  en¬ 
ables  one  to  isolate  subtle  effects  that  arise  specifically 
from  the  differences  between  the  two  systems.  Recent 
measurements  of  the  differential  cross  section  for  ion 
pair  formation  in  collisions  between  halogen  molecules 
and  alkali  atoms  or  metastable  rare  gas  atoms  provide 
such  an  opportunity.  Over  the  energy  range  30-150  eV, 
all  of  the  major  features  observed1  in  the  angular  dis¬ 
tribution  for  the  reaction 

Ar*(sPj,0)  +  Ij-Ar*+Ii  (1) 

can  be  correlated  with  those  observed2-*  in  the  analogous 
reaction 

K  +  Ij-K'-t-T,  .  (2) 

However,  there  are  differences  in  the  ratios  of  the  in¬ 
tensities  of  various  features  and  in  the  position  of  the 
rainbow  angle.  The  general  similarity  is  consistent  with 
the  fact  that  both  K  and  Ar*(2P)  have  an  active,  unpaired 
4s  electron.  The  observed  differences  between  the  two 
reactions  should  be  related  to  the  inner  shell  vacancy  in 
Ar*  and  to  the  additional  competing  channels  energeti¬ 
cally  accessible  to  the  system,  including  Penning  ioniza¬ 
tion: 

Ar*  +  l,-Ar  +  I*,  +  e-  .  (3) 

In  this  paper  we  present  model  calculations  designed  to 
elucidate  the  causes  of  the  observed  differences  in  the 
scattering  data. 

Ion  pair  formation  in  the  alkali -halogen  system  has 
been  extensively  investigated*- r  both  experimentally  and 
theoretically.  These  reactions  have  been  modeled  quite 
successfully  by  Los  and  co-workers*’*  using  classical 
trajectory  surface-hopping  techniques1'*  that  demon¬ 
strate  the  essential  features  of  the  reaction  mechanism 
leading  to  ion  pair  production.  Here  we  have  modified 
and  extended  those  model  calculations  to  Include  cou¬ 
pling  to  the  continuum  (Reaction  (3)1.  We  estimate  the 
autoionization  rate  rtf?)  using  a  simplified  form  of  a 
technique  proposed  recently  by  Miller  and  Morgner.  ** 


We  are  therefore  able  to  perform  quantitative  calcula¬ 
tions  that  examine  the  effect  of  competing  channels  upon 
the  ion  pair  distributions  in  Reaction  (1). 

In  Sec.  II  we  present  the  details  of  the  model  calcu¬ 
lations  and  discuss  the  inclusion  of  continuum  coupling. 
Our  results  are  presented  and  discussed  in  Sec.  IH, 
and  Sec.  IV  contains  a  summary. 

II.  THEORY 

A.  Summary  of  the  model 

We  have  adopted  the  theoretical  approach  developed 
by  Los  and  co-workers, 2’*’*  which  is  a  classical- trajec¬ 
tory,  surface-hopping  model*’*  with  the  following  addi¬ 
tional  approximations: 

1.  Simple  analytic  potentials  are  constructed  from 
pairwise  interactions  for  the  covalent  electronic  state 
K  + 12  and  for  the  ionic  state  K*  + 1 J. 

2.  The  cross  sections  are  obtained  as  the  average  of 
cross  sections  for  various  orientations  of  the  target 
molecule.  The  axis  of  this  target  is  assumed  not  to  ro¬ 
tate  during  the  collision.  Molecular  vibration,  however, 
is  included,  as  explained  below. 

3.  For  each  trajectory,  the  classical  vibration  of  the 
molecule  is  calculated  numerically  assuming  Morse  po¬ 
tentials,  and  the  deflection  of  the  projectile  is  obtained 
analytically  using  classical  perturbation  theory. 

An  essential  feature  of  the  model  is  that  the  ionic  and 
covalent  surfaces  cross  at  an  Ar*-It  distance  R  that  de¬ 
pends  on  the  vibrational  coordinate  of  the  molecule. 

This  fact  influences  the  collision  in  the  following  way. 

The  system  begins  asymptotically  on  the  covalent  sur¬ 
face.  When  the  particles  reach  the  location  of  the 
crossing  with  the  ionic  surface,  the  Landau-Zener 
formula11  is  used  to  compute  the  probability  of  switching 
to  the  other  surface.  As  shown  in  Fig.  1,  there  are  two 
possible  paths  to  ion  pair  formation.  The  surface  cross¬ 
ing  may  occur  on  the  inward  part  of  the  trajectory  (the 
dashed  line)  or  on  the  outward  part  (shown  by  the  dotted 
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FIG.  1.  Two  possible  paths  a  and  b  to  ion  pair  production  are 
schematically  indicated.  The  potential  curves  represent  cuts 
o(  the  Ar*  +  Ij  and  Ar*+  Ij  surfaces  at  a  fixed  I;  separation. 

line).  A  transfer  to  the  ionic  surface  will  initiate  an  ex¬ 
pansion  of  the  I j  bond  since  Ij  has  a  larger  equilibrium 
internuclear  distance  than  I}.  If  this  transfer  occurs  on 
the  incoming  trajectory,  the  resulting  I  j  expansion  will 
yield  an  effective  vertical  electron  affinity  that  increases 
with  time.  This  will  increase  the  radius  R  at  which  the 
second  surface  crossing  is  encountered  on  the  outward 
part  of  the  trajectory.  Since  the  matrix  element  con¬ 
necting  the  two  surfaces  decreases  strongly  with  in¬ 
creasing  R,  the  probability  of  a  diabatic  path  through  the 
second  crossing  may  be  significantly  larger  than  at  the 
first  crossing.  These  effects  have  already  been 
thoroughly  discussed  in  the  literature. 1-6 

The  analytic  form  of  the  potential  surface  we  have  used 
differs  somewhat  from  the  form  used  in  Ref.  2.  Ours  is 
defined  in  terms  of  Fig.  2  to  be 

F0#=F(R1)  +  F(RJ)  +  l.,(r)  ,  (4) 

Fu=F<R,)-^  +  F(Rj)-^5+u,(r).  (5) 

H  j 

[>o (r)  and  ot(r)  are  Morse  potentials  for  Ij  and  Ij,  respec¬ 
tively.  The  same  Morse  parameters  were  used  as  in 
Ref.  2,  namely,  0=1.54  eV,  0=1.87  A'1,  andr,  =  2.  67 
A  for  I,,  and  D'  =  1. 02  eV ,  0'  =  1.  23  A-1,  and  r',  =  3.  20  A 
for  Ij.  V{Rt)  is  a  repulsive  term  of  the  Born-Mayer 
form,  *’ 

V'«()=Ae-,“  i,  (6) 

where,  following  Ref.  2,  we  initially  use  A  =  8.  4  x  104 
eV  and  0  =  4.  762  A*1.  Note  that  for  simplicity  the  same 
constants  A  and  0  are  used  here  for  the  ionic  and  co¬ 
valent  terms.  The  potentials  we  use  are  much  simpler 
than  those  of  Ref.  2.  Since  ab  initio  calculations  are 
not  available,  we  choose  to  use  potentials  with  as  lew 
adjustable  parameters  as  possible. 

The  potential  surfaces  defined  by  Eqs.  (4)  and  (5)  do 
not  reflect  the  multiplicity  of  states  arising  from  the 
inner  shell  vacancy  in  Ar*.  Our  simplified  model  is 
equivalent  to  assuming  that  all  of  the  covalent  surfaces 
are  degenerate,  and  that  all  of  the  ionic  surfaces  are 
degenerate.  This  assumption  appears  reasonable  at 
large  values  of  R,  where  the  surfaces  cross.  In  this 
region  we  expect  no  inner  shell  rearrangement  pro¬ 
cesses.  The  core  electrons  should  have  only  a  minor 


influence  on  tne  potential  shapes  and  coupling  matrix 
elements.  This  fact  is  indeed  the  major  reason  for  the 
strong  analogy  to  alkali  reactions  and  the  justification 
for  our  use  of  the  single  surface  crossing  and  the  cou¬ 
pling  matrix  element  angular  dependence1’  [cosa  in  Eq. 

(7)  below]  that  are  appropriate  to  alkali  systems.  The 
inner  shell  vacancy  in  Ar*  would,  however,  be  expected 
to  influence  surface  parameters  at  much  smaller  R  val¬ 
ues  and  thereby  could  affect  the  shape  of  the  differential 
cross  sections.  We  will  return  to  this  point  later. 

We  use  the  same  form  for  the  coupling  matrix  ele¬ 
ment  as  in  Ref.  2: 

Ft2(R,r,9)=cI[/E.A.(r)]1/’R*e*'J**co8o  ,  (7) 

where 

R*=  (R/>/?)[f 1/4  +  E.  A.  (r)l/I1.  (8) 

I  is  the  ionization  potential  of  the  projectile  and 
E.  A.  (r)  is  the  vertical  electron  affinity  of  the  molecule 
as  a  function  of  the  separation  r.  Note  that  the  diagonal 
terms  FMand  V,,  are  orientation  dependent  because  they  are 
defined  in  terms  of  Ri  and  Rt.  Note  also  that  the  angular 
dependence  of  Vj2  forces  it  to  be  zero  at  the  angles  re¬ 
quired  by  symmetry. 15  Initially  the  constants  c(  =  0. 28 
a.u.  and  cj  =  0. 65  a.u.  were  used’  for  both  Reaction  (1) 
and  Reaction  (2). 

For  each  orientation  of  the  target  (in  a  space-fixed 
frame),  we  wish  to  calculate  the  two  branches  of  the 
classical  deflection  function.  These  two  branches  give 
the  scattering  angle  as  a  function  of  impact  parameter 
for  trajectories  in  which  the  surface  crossing  occurs  on 
the  inward  or  outward  part  of  the  trajectory.  (If  the 
crossing  occurs  twice,  or  not  at  all,  ion  pair  formation 
does  not  occur  and  the  trajectory  is  not  counted. )  For 
each  impact  parameter,  the  scattering  angles  are  there¬ 
fore  calculated  for  two  sequences  of  events.  In  both 
cases,  the  Ij  molecule  begins  at  rest  on  the  neutral  po¬ 
tential  curve  t>»(r)  at  the  equilibrium  distance  r„  and  the 
projectile  follows  a  straight-line  path  with  constant  ve¬ 
locity  throughout  the  collision.  For  the  ionic  branch  of 
the  deflection  function,  the  system  is  assumed  to  switch 
to  the  ionic  potential  at  the  first  crossing.  Then  the  Ij 
molecule  begins  to  move  classically  on  the  potential 
V\(r),  beginning  at  rest  at  r  =  r,  (I2).  The  location  of  the 
second  crossing  is  determined  numerically  by  monitor¬ 
ing  the  oscillatory  motion  of  the  target  and  the  rectilin¬ 
ear  motion  of  the  projectile.  For  the  covalent  branch 
of  the  deflection  function,  the  molecule  does  not  begin  to 
vibrate  until  the  crossing  is  encountered  a  second  time. 


I 

FIG.  2.  Illustration  of  the  coordinate  system  used  to  define 
the  potential  functions. 
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so  that  the  first  and  second  crossings  occur  at  the  same 
intermolecular  separation  R.  For  both  branches  of  the 
deflection  function,  the  orientation  angle  a  is  computed 
at  both  crossings.  The  orientation  influences  the  tran¬ 
sition  probability  through  the  factor  cosa  in  Eq.  (7). 

The  preceding  discussion  shows  how  the  appropriate 
potential  is  determined  at  each  point  along  a  given  tra¬ 
jectory.  The  scattering  angle  is  then  determined  by  in¬ 
tegrating  the  eumulative  momentum  transfer  perpendicu¬ 
lar  to  the  rectilinear  motion,  according  to  standard 
formulas'*  of  classical  perturbation  theory.  The  deflec¬ 
tion  is  easily  calculated  analytically  because  the  func¬ 
tions  in  Eqs.  (4)  and  (5)  are  written  as  sums  of  two- 
body  potentials. 

In  order  to  obtain  the  cross  section,  it  is  also  neces¬ 
sary  to  calculate  the  probability  that  the  system  will  fol¬ 
low  a  given  trajectory.  This  is  done  using  the  Landau- 
Zener  formula. 15  The  total  probability  for  a  given  se¬ 
quence  of  surface  crossings  or  avoided  crossings  is  the 
product  of  the  probabilities  for  each  event  separately. 
This  is  easily  calculated  once  the  locations  of  the  two 
crossing  points  are  known.  The  result  for  the  cross 
section  (for  a  particular  orientation  of  the  molecule)  can 
be  summarized  in  the  following  formula: 


/f{}'  is  the  matrix  element  for  the  one-electron  process 

Ar’  +  Ij-Ar  +  IJ  .  (14) 

Miller  and  Morgner10  pointed  out  that  H't\'  can  be  esti¬ 
mated  using  the  semiempirical  correlation  formula  of 
Olson  et  at. ": 

*GReR,c  ,  (15) 

where  R  is  the  internuclear  separation  in  a.u. ,  and 

1/c=vT(/‘* +/*»)  ’  (16) 

G=UhUr)ln/C  ,  (17) 

where  ft2  and  fAr  are  ihe  ionization  potentials  of  I2  and 
Ar  in  a.  u. ,  respectively.  I Htl  I  is  then  given  also  in 
a.  u. 

Based  on  past  experience1’  calculating  continuum  over¬ 
lap  factors  Sn,  we  have  simplified  the  Miller-Morgner 
formula  by  setting  S(i » 1.  This  approximation  enables 
us  to  evaluate  Htl  without  numerical  integration.  Fur¬ 
thermore,  since  the  potential  surfaces  are  not  in  gen¬ 
eral  known,  we  calculate  the  density  of  states  2xp  from 
the  energy  of  the  ejected  electron  as  R  -  «.  This  infor¬ 
mation  is  available  from  spectroscopic  data. 


7kain6=L 


2nbi  Pla'{bi) 
\d8l°'/db\ 


2nb,‘  Pcc,(bJ-) 


pio«(pco»)  is  u,e  overan  probability  of  following  an  ionic 
(covalent)  trajectory  for  a  particular  impact  parameter 
b.  The  sums  are  understood  to  be  over  all  impact  pa¬ 
rameters  that  lead  to  the  scattering  angles  8  or  -  8. 
Although  Eq.  (9)  presents  the  formal  definition  of  the 
cross  section,  we  found  it  useful  in  practice  to  use  a 
conventional  histogram  procedure. 


B.  The  Penning  ionization  channel 

1.  Estimating  T (R)  for  molecular  systems 

V{R)/K  is  the  rate  at  which  particles  are  lost  from  an 
initial  channel  into  the  ionization  channel  via  the  reac¬ 
tion  (3).  T  is  given  by"’1’ 

r=2*p|/f„|l  (10) 

=  2irp  \R |  |*  ,  (11) 

where  H  is  the  Hamiltonian,  *l<iHU  is  the  initial  bound 
state  wave  function  (Ar*  + 1?),  is  the  continuum 

configuration  (Ar  +  Ij  +  e‘),  and  p  is  the  density  of  con¬ 
tinuum  states.  If  the  continuum  wave  function  is  nor¬ 
malized  to  unit  amplitude  asymptotically,  then 


Miller  and  Morgner10  estimate  that  their  formula  is 
reliable  within  a  factor  of  3-5.  We  expect  that  our  addi¬ 
tional  approximations  will  degrade  the  accuracy  some¬ 
what,  but  the  result  should  certainly  still  be  a  reason¬ 
able  order -of -magnitude  estimate. 

2.  Ionization  probability  as  a  function  of  impact 
parameter 

The  autoionization  rate  obtained  in  the  previous  sec¬ 
tion  has  the  form  r(fi)oc  R1  e~M .  However,  we  have 
found  that  over  the  important  range  of  R  probed  in  the 
collisions,  the  numerical  values  of  T  can  be  well  fit 
(-10%)  by  the  simpler  form 

r(/?)  =  ot  .  (18) 

This  is  the  form  of  T  that  has  been  assumed  in  many 
previous  semiempirical  studies.  20,21 

We  have  found  that  for  a  T  of  the  form  of  Eq.  (18),  a 
very  simple  formula  can  be  obtained,  in  the  perturba¬ 
tion  limit,  for  the  ionization  probability  as  a  function  of 
impact  parameter.  Since  the  translational  motion  is  al¬ 
ready  treated  using  classical  perturbation  theory,  a 
perturbation  assumption  for  the  continuum  coupling  is 
consistent  with  the  rest  of  the  model. 


2irp  =  4/fe  ,  (12) 

where  the  kinetic  energy  of  the  ejected  electron  is,  in 
atomic  units,  £  =  1/2  fe2. 

Miller  and  Morgner10  have  suggested  that  the  matrix 
element  Htt  can  be  written  as  the  product  of  an  overlap 
factor  and  a  charge  exchange  factor.  Soecializlng  to 
our  case,  Ar*  +  Ij,  this  would  be  the  product 

/f(larSIlxff{I*  (13) 

where  S„  is  the  overlap  between  the  valence  electron 
orbital  on  Ar*  and  the  final  state  continuum  orbital,  and 


Let  6  be  the  impact  parameter.  Then  the  survival 
probability  P(b)  that  ionization  does  not  occur  on  a  tra¬ 
jectory  of  impact  parameter  b  is 

P(b)=  exp(-  2F(6)|  ,  (19) 

where  in  the  perturbation  limit 


T  (R)dR 


(20) 


and  vt  is  the  incident  velocity. 

Substituting  Eq.  (18)  into  Eq.  (20),  integrating  by 
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parts,  and  substituting  x  =  R/b,  we  obtain 

F(b)=  f~G*^ I  e‘*'dx  .  (21) 

The  integral  can  be  represented  in  terms  of  the  modi¬ 
fied  Bessel  function  Kt.  Then 

F(6)=  (atb)  .  (22) 

Further  simplification  is  possible  because  of  the  relation 

xKx  Cv)*2.076e°-,m  (23) 

which  is  valid  to  about  10%  for  2s.vs  12.  The  final  re¬ 
sult  is 

(24) 

ajftvn 

Formulas  (19)  and  (24)  may  now  be  incorporated  into 
the  model  in  the  following  way.  For  each  orientation,  a 
deflection  function  is  calculated  that  gives  the  angular 
scattering  as  a  function  of  impact  parameter.  When  the 
cross  section  is  calculated  from  a  particular  branch  of 
the  deflection  function,  a  factor  is  first  included  to  rep¬ 
resent  the  probability  that  the  particles  follow  the  cor¬ 
rect  sequence  of  curve  crossings.  The  result  is  then 
multiplied  by  P(b),  which  gives  the  probability  that  the 
trajectory  is  completed  without  loss  to  ionization. 

III.  RESULTS  AND  DISCUSSION 

A.  K+lj 

This  reaction  has  already  been  treated  in  detail. 2-4 
Here  we  demonstrate  that  the  use  of  the  potentials  de¬ 
fined  in  Sec.  HA  leads  also  to  satisfactory  agreement 
with  experiment.  In  Fig.  3,  we  compare  the  calculated 
and  experimental  values  at  a  lab  energy  (£a)  of  60  eV. 

The  agreement  is  not  surprising  since  our  potentials  are 
nearly  the  same  as  those  fitted  to  the  data  in  Ref.  2. 

The  major  difference  is  that  we  have  neglected  the  small 
van  der  Waals  term. 

B.  Ar*  +  I2:  General  approach 

Gillen  el  al. 1  pointed  out  that  their  results  for  Ar*-+I2 
were  similar  to  those  of  Aten  et  al. 2  for  K  +  I2.  The  im¬ 
portant  differences  are  that  the  covalent  peak  is  rela¬ 
tively  somewhat  smaller  for  Ar*  + I2  than  for  K  + 12  and 
that  the  rainbow  in  the  ionic  peak  occurs  at  smaller 
angles.  The  first  difference  could  be  caused  :y  the  con¬ 
tinuum  coupling  in  the  metastable  system  or  by  differ¬ 
ences  in  the  potentials  and  charge  exchange  matrix  ele¬ 
ments.  The  second  difference  should  be  associated  with 
a  difference  in  the  shape  of  the  potential  surfaces. 

In  this  and  the  following  subsections  we  describe  our 
quantitative  calculations  to  investigate  possible  causes 
for  the  experimental  differences.  With  only  a  qualita¬ 
tive  knowledge  of  the  Ar*  + 12  potential  surfaces,  we  do 
not  expect  to  determine  uniquely  the  coupling  parameters 
in  this  system.  Hence,  our  approach  to  Ar*  +I2  ion  pair 
formation  calculations  is  to  start  with  the  K  +  I2  parame¬ 
ters  and  determine  the  changes  in  the  distributions 
caused  by  individually  altering  the  continuum  width  (r 
=  0  initially),  the  surface  coupling  "arameters,  and  the 


r  leV-degrees) 


FIG.  3.  Comparison  of  previous  calculations  and  experimental 
data  (Ref.  2)  on  the  system  K  +  U  with  the  present  calculation, 
which  uses  slightly  different  interaction  potentials. 


interaction  potentials.  The  results  should  yield  insight 
into  the  possible  causes  of  the  differences  observed  be¬ 
tween  Ar*  + 12  and  K  + 12  in  the  ion  pair  channel. 

We  initially  performed  trial  calculations  for  Ar*  +  I2 
assuming  r  =  0,  and  using  the  same  potential  parameters 
as  for  K  +  I2.  The  crossing  point  differs  slightly  for  this 
system  because  the  ionization  potential  of  Ar*  is  slightly 
smaller  than  that  for  K.  This  change  increases  Rc 
somewhat,  and  thereby  decreases  V,2  slightly,  but  the 
effect  on  the  calculations  was  quite  small.  The  calcu¬ 
lated  peak  ratios  and  the  position  of  the  rainbow  angle 
were  essentially  the  same  as  those  calculated  for  K  +  I2, 
and  these  values  did  not  agree  with  the  data  for  Ar*  + 12. 

.  C.  Ar*  +  I2:  The  continuum  channel 

We  then  examined  the  effect  of  including  continuum 
coupling.  It  should  be  noted  that  T(R)  is  assumed  to  af¬ 
fect  only  the  covalent  potential.  This  is  consistent  with 
the  conclusions  of  Hultsch  et  al. 22  on  similar  systems. 
One  argues  that  Penning  ionization  of  the  covalent  elec¬ 
tronic  state  is  initiated  by  a  one-electron  process, 
namely,  the  "exchange”  process,  in  which  an  electron 
from  the  ground  state  atom  jumps  into  the  hole  of  the 
excited  atom.  There  is  then  a  large  overlap  between  the 
valence  electron  of  the  excited  atom  and  the  continuum 
orbital.  In  contrast  for  a  system  on  the  Ar*  +I2  surface, 
continuum  coupling  necessitates  a  simultaneous  two- 
electron  rearrangement;  such  a  process  would  be  ex¬ 
pected  to  have  a  much  lower  probability. 

For  the  ground  electron  state  of  I2,  which  is  9. 3  eV 
above  I2,  we  foll'w  the  prescription  of  Sec.  n  B  and  ob¬ 
tain  (for  T  and  R  in  atomic  units) 

r(R)^1.752R,e-,-HM  .  (25) 
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FIG.  4.  Calculation  of  the  differential  croea  section  for  ion 
pair  formation  in  Ar*  *■  I,  compared  with  the  data  of  Ref.  1. 

The  effect  of  different  autolonizatlon  rates  T(R)  In  the  covalent 
channel  is  Illustrated.  r0  Is  the  function  estimated  using  the 
method  discussed  in  Sec.  n.  B.l,  and  given  by  Eq.  (261.  Suc¬ 
cessively  larger  values  of  HR)  reduce  the  covalent  peak  by 
larger  fractions. 


For  R  in  the  range  4  -  12a„,  this  is  well  fit  by  the  form 
r„«)arll.  16e*MU"  a.u.  (26) 

We  performed  calculations  using  this  form  of  r(ft)  and 
keeping  the  other  potential  parameters  the  same.  The 
result  is  that  the  continuum  coupling  is  not  large  enough 
to  cause  a  significant  effect  on  the  ion  pair  distributions. 

The  charge  transfer  matrix  element  estimated  from 
Eq.  (15)  is  a  one-electron  matrix  element.  In  fact,  r(R) 
should  be  larger  than  r0W)  to  account  for  the  two  equiv¬ 
alent  electrons  in  the  2ir} n<,  orbital.  Moreover,  the  I2 
molecule  has  a  valence  configuration  (a, 5/>)2  (itin.ufy)1 
(’rJ/s,.5P)2  (*« and  these  ten  electrons 
have  vertical  ionization  potentials2’  between  9.  3  and 
13. 0  eV ,  In  the  spirit  of  our  one-electron  estimate  of 
r„,  we  should  also  calculate  approximate  continuum 
coupling  widths  for  the  other  electrons. 

Although  the  ionization  potential  for  electrons  In  the 
lowest  two  of  these  orbitals  Is  larger  than  the  Ar*  ex¬ 
citation  energy,  Miller  and  Morgner10  have  shown  that 
the  electronic  energy  transfer  to  a  pseudocontinuum 
Rydberg  state, 

Ar*  +  I2-Ar  +  If*  (27) 

(core  excited  In  this  case)  can  be  treated  in  an  equivalent 
way  to  a  true  continuum  interaction.  The  details  of  the 
formula  are  slightly  different;  however,  the  calculated 
coupling  widths  are  found  to  be  quite  similar  to  those 


for  states  whose  ionization  potentials  are  lower  than  the 
excitation  energy  of  the  incoming  state.  Assuming  the 
contributions  for  all  ten  electrons  can  be  added  (a  rea¬ 
sonable  first  approximation),  we  obtain  a  resultant  r (ft) 
that  is  10-15  times  larger  than  r9(ft)  for  R  values  be¬ 
tween  5  and  11  a.u. 

The  effects  of  various  coupling  functions  T(R)  on  the 
calculated  angular  distribution  for  Reaction  (1)  at  28  eV 
are  shown  in  Fig.  4.  Continuum  effects  are  generally 
more  pronounced  at  lower  energies,  and  Et  =  28  eV  is 
the  lowest  energy  studied.  Clearly  T=  10  r,,  which  ap¬ 
proximates  the  summed  ten  electron  continuum  coupling 
extimated  above,  yields  a  noticeable  depletion  of  the  co¬ 
valent  peak.  However,  r>lOOr0  seems  required  before 
the  intensity  ratios  of  the  covalent  and  ionic  peaks  are  in 
reasonable  agreement  with  experiment.  Since  our  ap¬ 
proximate  estimates  of  T  are  relatively  crude,  it  is  pos¬ 
sible  that  r  is  indeed  an  order  of  magnitude  larger  than 
estimated.  Hence,  the  continuum  coupling  effect  could 
make  some  contribution  to  the  attenuation  of  the  covalent 
peak  in  Reaction  (1).  However,  the  larger  T  values 
(100  I'd,  1000  r„)  give  covalent  peaks  that  are  strongly 
attenuated  at  low  r  values;  the  angular  positions  of  the 
calculated  peaks  then  disagree  considerably  with  the 
data. 

O.  Ar*  + 12:  Sensitivity  to  the  coupling  matrix  element 

Another  possible  cause  for  the  smaller  covalent  peak 
in  Ar*  +  I2  could  be  a  larger  coupling  matrix  element  V,2. 
We  have  performed  calculations  in  which  T  =  0  and  the 
constant  c,  of  VXJ  is  varied  [see  Eq.  (7)].  We  find  that  a 
significant  change  in  from  the  value  for  K  +  I2  is 
needed  to  cause  the  desired  change  in  Ar*  +I2  peak 
ratios.  Figure  5  shows  the  experimental  angular  distri¬ 
bution  at  £g  =  51  eV  ,  compared  with  three  calculations 
assuming  K,2  multiplied  by  factors  of  1.  0,  2.  5,  and  5.0 
relative  to  K  +  I2.  Even  though  the  peak  ratios  are  accu¬ 
rately  matched  by  scaling  F,2  by  approximately  2.  5,  we 
feel  that  this  is  an  unreasonable  change  considering  the 
expected  similarities  between  K  and  Ar*.  Although  the 
potential  surfaces  for  the  two  reactions  have  quite  dif¬ 
ferent  symmetries,  these  differences  are  associated 
with  the  symmetry  of  the  Ar*  core.  At  the  very  large 
crossing  radius  (-5.7  A)  between  the  ionic  and  covalent 
surfaces,  the  properties  of  the  core  should  play  a  minor 
role  in  the  electron  transfer  process  and  the  outer  Ar* 
electron  should  behave  like  a  4s  electron  in  K. 

E.  Ar*  +  I2:  Sensitivity  to  V'22 

It  is  obvious  from  Fig.  5  that  the  calculated  rainbow 
angle  is  not  in  good  agreement  with  the  experimental 
data.  None  of  the  modifications  detailed  above  has  al¬ 
tered  the  calculated  rainbow  position.  Since  its  position 
is  most  sensitive  to  the  well  depth  in  the  ionic  channel, 
we  Investigated  modifications  to  the  potential  parameters 
necessary  for  achieving  a  match  to  the  experimental 
rainbow  angle.  In  our  model,  the  well  depth  Is  con¬ 
trolled  by  the  choice  of  the  parameters  A  and  B  in  Eq. 

(6).  We  found  that  varying  these  parameters  to  decrease 
the  well  depth  caused  the  rainbow  to  shift  to  smaller 
angles,  as  expected.  Another  effect  of  this  change  was 
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FIG.  5.  Calculation  of  the  differential  crosa  section  for  ion 
pair  formation  in  Ar*  + 1:  compared  with  the  data  of  Ref.  1. 
The  effect  of  different  coupling  matrix  elements  is  Illus¬ 
trated.  run  Is  set  to  zero  in  all  cases.  la  the  form 
given  by  Eq.  <7>  with  c(»0.28  and  c2>0.65,  the  same  param¬ 
eters  used  In  Ref.  2  for  K  ♦  f,.  The  scaled  values  of  F,2  cor¬ 
respond  to  changing  c,  to  0.7  and  1.4. 


a  reduction  in  the  size  of  the  covalent  peak  relative  to 
the  ionic  peak.  Figure  6  shows  the  effect  of  changes  in 
the  potential  parameters  upon  the  ion  pair  distributions. 
The  modification  that  gave  a  reasonable  fit  to  the  ionic 
rainbow  angle  corresponds  to  an  Ar*-IJ  “well  depth” 
that  is  -0.  5  eV  shallower  than  that  of  K  *  + 1  j  (at  the 
same  I-I  internuclear  separation). 

For  comparison,  there  are  good  ab  initio  calculations 
for  the  diatomic  potentials  in  ArF,  KrF,  and  the  Xe 
halides. 24  For  each  of  these  systems  the  well  depths  of 
the  various  ion  pair  states  are  between  0. 3  and  0. 7  eV 
shallower  than  those  of  the  analogous  alkali  halide  mole¬ 
cule.  With  spin-orbit  cot  pling  included,  the  three  ion 
pair  states  calculated  for  each  rare-gas  halide  pair 
differ  by  less  than  -0. 2  eV  in  well  depth  and  0. 1  A  In 
R,  from  each  other.  If  a  similar  result  obtained  for  the 
trlatomic  systems,  our  method  of  treating  the  Ar*  +  1, 
system  by  a  single  (average)  Ion  pair  surface  intersect¬ 
ing  a  single  covalent  surface  would  be  a  reasonable  ap¬ 
proximation  to  the  average  scattering  on  the  relevant 
surfaces  and  our  resultant  0.  5  eV  difference  in  Ar*-Ij 
and  K’-I  j  well  depths  would  appear  to  be  reasonable. 

F.  Ar*  + 12:  Conclusions 

Our  calculations  clearly  indicate  that  the  ion  pair  an¬ 
gular  distributions  for  Reaction  (1)  depend  quite  sensi¬ 


tively  on  the  shapes  of  the  Interaction  potentials,  the 
surface  coupling  matrix  elements,  and  the  amount  of 
depletion  from  competing  continuum  coupling  and  exci¬ 
tation  transfer  channels.  With  only  crude  model  poten¬ 
tials  and  order -of -magnitude  estimates  for  the  continuum 
coupling  and  excitation  transfer  effects,  it  does  not  seem 
Justified  to  vary  parameters  to  achieve  a  fit  to  the  data 
that  would  be  both  model  dependent  and  nonunique.  How¬ 
ever,  an  estimate  of  the  probable  causes  for  the  differ¬ 
ences  from  Reaction  (2)  can  be  useful  in  suggesting  di¬ 
rections  for  further  experimental  or  theoretical  effort 
on  this  system. 

The  shift  of  the  ionic  rainbow  for  Reaction  (1),  to 
lower  r  values  than  that  for  Reaction  (2),  seems  to  be 
consistent  with  a  shallower  well  in  the  Ar*  +1}  interac¬ 
tion  (see  Fig.  6).  This  change  also  lowers  the  magni¬ 
tude  of  the  covalent  peak  relative  to  the  ionic  peak.  In 
fact  all  other  modifications  combined  would  only  need  to 
lower  the  covalent  peak  by  an  additional  30%-4O%  to 
achieve  a  good  agreement  with  the  experimental  data  at 
£0=51  eV.25 

Combined  with  the  change  in  the  ionic  potential  de¬ 
scribed  in  Sec.  E,  an  increase  of  c(  by  -40%  over  the 
K  + 1,  value  (0. 28)  would  achieve  a  reasonable  match 
to  the  experimental  data.  However,-  such  a  change  in 
F,2  still  seems  rather  large.  As  evidence  for  this  as¬ 
sertion,  we  note  that  Hubers  et  at. 24  present  a  reduced 
matrix  element  formula  that  is  slightly  modified  from 
Eq.  (7)  and  correlates  with  all  alkali-halogen  systems  to 
an  accuracy  of  better  than  10%.  The  modification  of  the 


FIG.  6.  Calculation  of  the  differential  croee  section  for  ion 
pair  formation  in  Ar*  «• I,  compared  with  the  data  of  Ref.  1. 
The  potential  parameters,  defined  in  Eqs.  (4)-(6),  are  A«8.4 
x  10*  eV  and  0 •4.762  A*1  for  the  original  well,  and  A  -1.35 
xio4  eV  and  0-3.781  A*1  for  the  shallower  well. 
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Ar*  +  Ij  matrix  element  postulated  here  would  deviate  by 
40%  from  the  prediction2'  based  on  the  alkali  data  and 
would  then  suggest  some  real  differences  of  Ar*  from 
alkalilike  properties  at  distances  of  R„  *  5. 7  A.  It  is 
difficult  to  understand  such  a  large  effect  at  that  distance 
caused  by  a  vacancy  in  the  ion  core. 

Alternatively,  a  continuum  coupling  width  of  ~  200  r0 
could  be  combined  with  the  shallower  ionic  potential  to 
match  the  peak  ratios  adequately.  This  width  is  some¬ 
what  larger  than  the  rough  estimate  (10-15  r,)  based  on 
the  Miller -Morgner 12  arguments,  and  the  attenuated  co¬ 
valent  peak  would  be  located  at  larger  r  values  than  the 
data  would  allow.  Yet  our  estimate  of  r  is  subject  to 
considerable  uncertainty  and  could  be  significantly  in 
error.  Also,  the  covalent  peak  position  and  shape  are 
sensitively  dependent  on  the  potential  functions  used. 

This  has  been  demonstrated4  in  K  +  Br2  by  a  comparison 
of  calculations  on  diabatic  and  adiabatic  surfaces. 

Hence,  continuum  coupling  (3)  and  excitation  transfer 
(27)  reactions  involving  the  ten  valence  Ij  electrons 
could  be  responsible  for  all  or  part  of  the  remaining 
30%-40%  depletion  of  the  covalent  peak. 

Finally,  we  note  that  making  the  potentials  more  re¬ 
pulsive  tends  to  decrease  the  size  of  the  covalent  peak 
relative  to  the  ionic  peak.  For  simplicity,  we  have  con¬ 
strained  our  calculations  so  that  the  repulsive  Born- 
Mayer  terms  are  identical  on  the  ionic  and  covalent  sur¬ 
faces.  Although  this  is  a  reasonable  first  approxima¬ 
tion,  one  could  decrease  the  covalent  peak  by  relaxing 
this  constraint  and  allowing  more  repulsion  in  the  co¬ 
valent  potential.  If  such  a  difference  were  responsible 
for  the  smaller  covalent  peak  in  (1)  relative  to  (2),  it 
would  necessarily  imply  differences  between  the  two 
systems  in  the  repulsive  region.  This  is  easier  to  ac¬ 
cept  than  large  differences  (40%)  in  the  coupling  matrix 
elements  at  the  much  larger  distances  associated  with 
the  surface  crossings. 

IV.  SUMMARY 

The  small  differences  between  the  differential  cross 
sections  for  ion  pair  formation  in  the  Ar*  +  Ij  and  K  +  Ij 
system  can  be  explained,  in  part,  by  a  more  repulsive 
wall  in  the  ionic  Ar*  +  Ij  surface  relative  to  K*  + 1  j. 

This  yields  a  smaller  ionic  well  depth  (by  -0.  5  eV  at  r, 
of  Ij)  and  a  correspondingly  smaller  ionic  rainbow  angle 
that  is  in  agreement  with  the  data.  If  the  covalent  sur¬ 
faces  Ar*  +  Ij  and  K  +  Ij  have  repulsive  terms  similar  to 
those  on  their  ionic  counterparts,  the  calculated  ratio 
for  the  intensity  of  the  covalent  peak  to  the  ionic  peak  in 
Reaction  (1)  is  too  large  by  a  factor  of  -1.4. 

This  ratio  can  be  improved  by  one  or  more  of  the  fol¬ 
lowing:  a  relative  increase  in  the  repulsion  on  the  Ar* 

+  Ij  covalent  surface;  a  depletion  of  the  covalent  peak  by 
competing  Penning  (3)  and  excitation  transfer  (27)  reac¬ 
tions;  and  an  increase  in  the  coupling  matrix  element  at 
the  surface  crossing  between  the  ionic  and  covalent  sur¬ 
faces.  The  third  possibility  seems  least  likely  since  it 
requires  effects  from  the  Ar*  core  vacancy  to  be  impor¬ 
tant  at  large  R  (-  5. 7  A). 

Since  the  potential  surfaces  have  only  been  estimated 
in  a  crude  way  and  since  the  continuum  coupling  width  is 


only  known  very  approximately,  we  feel  that  no  claim  to 
a  unique  fit  can  be  made.  Continuum  coupling  and  exci¬ 
tation  transfer  processes  may  indeed  cause  observable, 
but  small,  depletions  of  the  covalent  peak  in  the  ion  pair 
process  (1). 
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An  analytic  approximation  to  the  differential  cross  section  is  derived  for  ion-pair  formation  in  a  model  two-state 
(“turned-on  Coulomb”)  system.  The  result  gives  insight  into  the  structure  observed  in  recent  experiments. 


1 .  Introduction 

Ion-pair  formation  differential  cross  sections  have 
been  measured  recently  for  fast  collisions  of  alkali 
atoms  [  1 1  and  metastable  rare  gas  atoms  [2]  with 
various  targets.  The  results  usually  have  similar  basic 
features  that  result  from  a  crossing  between  the 
incoming  covalent  potential  surface  and  a  strongly 
attractive  ion-pair  surface.  The  ion-pair  product 
angular  distribution  generally  exhibits  two  peaks,  the 
“ionic”  and  the  “covalent”  peaks,  resulting  from  tra¬ 
jectories  that,  at  internuclear  separations  smaller  than 
the  surface  crossing  distance,  travel  along  the  ionic 
and  covalent  surfaces,  respectively.  The  ionic  peak, 
located  at  larger  angles,  often  has  a  single  secondary 
maximum  or  shoulder  on  its  low-angle  side.  This 
secondary  structure  is  observed  with  many  molecular 
targets  [  1 ,2] ,  but  also  appears  in  the  simple  two-body 
reaction  Na  +  I  -»  Na*  +  1*  (3) .  In  this  paper  we  pre¬ 
sent  a  very  simple  model  that  highlights  the  origin  of 
this  secondary  maximum  and  suggests  why  it  appears 
in  some  systems  but  not  in  others. 

2.  Model 

In  this  section  we  derive  an  analytic  approximation 
to  the  differential  cross  section  for  ion-pair  formation 
under  idealized  assumptions  about  the  two  electronic 
potentials  involved.  We  find,  a  posteriori,  that  the 
trajectories  that  contribute  to  the  secondary  maximum 
only  probe  the  potentials  in  a  range  where  our  approx¬ 


imations  are  reasonably  accurate. 

Consider  the  model  potentials  shown  in  fig.  1 .  We 
assume  that  the  ionic  potential  is  given  by  the  pure 
Coulomb  form  F,  (R)  =  - 1/7?.  Let  A£  be  the  energy 
difference  between  the  neutral  (covalent)  and  ion-pair 
channels.  Then  if  the  zero  of  energy  is  defined  to  be 
the  asymptote  of  the  ionic  potential,  the  covalent 
curve  is  approximated  by  the  constant  value  F0(R)  = 


Fig.  1.  Model  potential!  and  curve  crowing  sequence  for  the 
time-reversed  ionic  trajectory  that  converts  ion-pair  reactants 
to  neutral  products. 
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-AE  *  -IIRC,  where  Re  is  the  crossing  point  of  K0(R) 
and  ('((ft).  (Atomic  units  are  used  throughout  this 
section  unless  otherwise  indicated.)  For  these  idealized 
potentials,  the  deflection  function  (scattering  angle 
versus  impact  parameter)  can  be  determined  analyti¬ 
cally  for  trajectories  that  begin  asymptotically  on  V0 
and  end  up  on  V,  after  a  curve  crossing  at  Rc.  The 
deflection  function  has  two  brancnes,  corresponding 
to  trajectories  that  make  the  crossing  on  the  inward 
or  outward  part  of  the  classical  motion. 

It  is  convenient  to  begin  by  considering  the  time- 
reversed  trajectory  shown  by  the  arrows  in  fig.  1 .  Let 
E  and  b  denote  the  cm.  collision  energy  and  impact 
parameter  relative  to  the  ionic  asymptote.  We  denote 
by  E'  and  b'  the  corresponding  quantities  relative  to 
the  covalent  asymptote.  Clearly  E1  ■  E  +  A E\  the  rela¬ 
tion  between  b  and  b'  is  obtained  by  invoking  the  fact 
that  the  angular  momentum  must  be  the  same  at  t  = 

+“>  and  r  *  -<*>.  The  result  is 

b’~b[\-(A^/E'Yn-  (1) 

Using  the  potentials  shown  in  fig.  1,  the  c.m. 
scattering  angle  0  for  an  impact  parameter  b  is  given 
by  (4) 


(dfl/R’Xl  ~  b'/R*  *  1  /«c£^~,/a . 


(2) 


The  turning  point  R0,  which  depends  on  b,  is  ob¬ 
tained  from  the  relation 

1  -*7*o  +  l/R0E*0,  (3) 

and  b<(R\  *  RJE)in  is  necessary  to  reach  Re . 
After  some  algebraic  and  trigonometric  manipulation, 
eq.  (2)  can  be  solved  exactly  for  the  deflection  func¬ 
tion 

0  *  -arcsin  Y  -  arcsin[Y(l  -  b2t/)l/2J ,  (4) 

whete 

y*(l+4b2£2r'/2  (5) 


and 

t/=(R2  +RJEfl.  (6) 

The  first  term  in  eq.  (4)  is  exactly  half  the  Rutherford 
deflection  angle  0r  for  pure  coulombic  scattering 

0r  *  -2  arcsin  Y  *  -2  arctan(l/2fl£).  (7) 

The  alternative  covalent  path,  which  involves  a  trans¬ 
fer  from  the  coulombic  to  the  flat  potential  on  the 
inward  portion  of  the  trajectory,  has  a  deflection 
angle  9'  given  by 

0'*0r-0,  (8) 

so  that  the  covalent  and  ionic  branches  of  the  deflec¬ 
tion  function  yield  scattering  angles  at  every  b  that 
add  up  to  Sr  for  the  same  b.  Considering  the  limiting 
case  where  the  two  branches  of  the  deflection  func¬ 
tion  meet  for  a  turning  point  R0(f>‘)  =  Rc,  one  ob¬ 
tains  from  eq.  (3) 

b*  »  U'xn 

and  the  second  term  in  eq.  (4)  vanishes,  implying  that 

0*  *  -arcsin  Y  =  \t)R  .  (9) 

Alternatively,  eq.  (8)  and  the  requirement  that  0’ ■  0 
yield  the  same  result.  Using  eq.  (3)  this  can  be 
reduced  to 

0*  *  -arcsin  [1/(1  +  2RCE)]  .  (10) 

One  can  show  that  this  result  is  in  agreement  with  a 
formula  given  by  Young  et  al.  [5)  for  the  same  limit¬ 
ing  angle. 

A  single  approximation,  the  small-angle  approxi¬ 
mation  tan  0  *  0,  is  all  that  is  required  to  invert  the 
deflection  function  and  obtain  the  differential  cross 
section.  Eq.  (4)  becomes 

0  *  -()/2MT)[l  +(1  — h1!/)1'2],  (11) 

and  the  differential  cross  section  for  the  inverse  of  ion 
pair  formation  (do0i/dw)  sin  0  is  given  by 

(do/dw)  sin  0  =  lb  db/d0| 

*|16£,20(f/-4£,202)/(4£,202  +£/)3l.  (12) 

Surprisingly  the  identical  eq.  (12)  is  obtained  for  the 
covalent  branch  of  the  deflection  function.  Note  that 
as  Rc  -»  the  limit  of  purely  coulombic  scattering, 
eq.  ( 12)  would  yield 

(do/dw) sin  0  -*  I/£20*,  (13) 
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which  agrees  with  the  Rutherford  result  under  the 
same  small  angle  approximation  6  »  tan  6  **  sin  6. 

Eqs.  (12)  and  (1)  can  be  used  to  obtain  the  differ¬ 
ential  cross  section  for  the  ion-pair  formation  process 
(which  will  be  called  “turned-on  Coulomb  scatter¬ 
ing”)  . 

(do,0/dcj)  sin  6  =  | b'  db'lde  |  (14) 

=  (1  -  Af/f1') [16£j0/(4£-20j  +  i/)33K/-  4£-10,|. 

Differentiation  shows  that  this  function  has  two 
peaks  at  reduced  deflection  angles 

r,  =E'd  =0.l81£'l/l'J/£  (IS) 

and 

rj  =0J9(£‘Ul,1IE .  (16) 

The  function  drops  to  zero  between  these  peaks  when 
\U  —  4£292l  =  0  and 

rmin=0  .StU^/E.  (17) 

Eq.  (17)  is  an  alternative  small  angle  approximation 
to  eq.  (10),  and  locates  the  minimum  in  the  differ¬ 
ential  cross  section  between  the  ionic  and  covalent 
branches.  The  peak  Ti  will  be  seen  to  approximate  the 
position  of  the  covalent  peak,  while  r2  gives  a  simple 
first  approximation  to  the  position  of  the  secondarv 
maximum.  Clearly,  the  ionic  rainbow  [1]  maximum 
cannot  be  obtained  in  our  simple  model  that  neglects 
repulsive  forces  entirely. 

By  solving  eq.  (1 1)  for  A,  one  can  demonstrate  that 
the  coulombic  impact  parameter  A2  that  corresponds 
to  the  peak  in  the  differential  cross  section  at  r2  is 

bt  *0.901  U'tn.  (18) 

Then  the  corresponding  impact  parameter  along  the 
flat  asymptote  is 

b'i  *  0.901Re.  (19) 

The  analogous  result  for  the  covalent  peak  is 
b\  =*0.64l/?c.  (20) 

3.  Example:  Na  +  I 

Delvigne  and  Los  [3]  have  measured  differential 
cross  sections  for  the  reaction 

Na  +  I  -*  Na*  ♦  I"  (21) 


Fig.  2.  Nat  potentials  of  Delvigne  and  Los  ( - )  compared 

to  our  model  potentials  ( - )  for  the  same  system. 


and  have  calculated  deflection  functions  for  this  sys¬ 
tem  us.i.g  a  multi-term  potential  function.  Their 
potentials  are  shown  in  fig.  2  along  with  idealized 
curves  consistent  with  our  model  developed  in  sec¬ 
tion  2.  The  major  differences  in  the  realistic  poten¬ 
tials  are,  of  course,  the  repulsive  terms  in  both  the 
ionic  and  covalent  curves,  but  there  are  R  ,  R , 
and  R'1  attractive  terms  in  their  Rittner  Na*  ♦  I" 
potential  that  cause  it  to  be  more  attractive  for 
RZ  3  A  and  move  the  curve  crossing  out  to  Rt  *  6.98 
A  compared  to  the  coulombic  Re  *  6.94  A  for  the 
simpler  approximation. 

Fig.  3  shows  deflection  functions  for  reaction  (21) 
at  £c  =  13.1  eV.  Results  of  Delvigne  and  Los  are 
compared  to  our  calculations  using  eqs.  (4)  and  (8). 
The  coulombic  case  [from  eq.  (7)]  is  also  shown  for 
comparison.  Here  our  results  deviate  from  those  of 
Delvigne  and  Los  in  ways  that  are  easily  associated 
with  the  differences  in  the  potentials  shown  in  fig.  2. 

The  measured  differential  cross  sections  for  reac¬ 
tion  (2 1 )  at  two  energies  are  shown  in  fig.  4  along 
with  our  curves  calculated  from  eq.  (14).  The  Ruther- 
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Fig.  3.  Deflection  functions  for  reaction  (21)  at  E'  ■  13.1  eV. 
The  results  of  Delvignc  and  Los  (•■  )  are  plotted  versus  the 

covalent  impact  parameter  b‘ .  Our  results  from  eqs.  (4)  and 
(8)  ( - )  have  impact  parameters  b'  and  b  at  large  dis¬ 

tances  along  the  flat  and  coulombic  potentials,  respectively. 
Shown  for  comparison  is  a  pure  coulombic  deflection  func¬ 
tion  plotted  versus  b  from  eq.  (7)  ( . )  for  —  «•, 

ford  cross  sections  derived  from  eq.  (13)  are  also 
shown  for  comparison.  Clearly  the  tumed-on  Coulomb 
calculation  cannot  produce  the  ionic  branch  rainbow 
at  t  w  240  eV  deg  since  repulsive  terms  are  not 
included.  Yet  the  calculation  at  38.7  eV  matches  all 
other  major  features  of  the  experimental  results,  in 
particular  the  shape  and  position  of  the  intermediate 
secondary  peak.  This  may  seem  surprising,  since  the 
calculation  does  not  take  into  account  the  probability 
P(b')  for  crossing  from  one  curve  to  another  at  the 
point  Rc.  This  neglect  is  in  effect  equivalent  to  an 
assumption  that  the  probability  of  following  either 
path,  ^(1  -  P)  or  (1  -  P)P,  is  nearly  constant  for  all  b' 


Fig.  4.  Smoothed  differential  cross  section  measurements  for 

reaction  (21)  from  Delvigne  and  Los  ( - )  compared  to 

our  turned-on  Coulomb  calculations  ( - )  and  purely 

coulombic  (Rutherford)  results  (•••—• — )  under  the  small 
angle  approximation  «  ■  sin  »  «  tan  «.  The  calculations  have 
not  been  smoothed  to  account  for  apparatus  resolution 
parameters.  The  abscissa  r  is  the  product  of  the  collision 
energy  and  the  scattering  angle. 

values  contributing  to  the  intermediate  peak  region  of 
the  differential  cross  section.  Coincidentally,  this 
seems  to  be  the  case  [3  J  at  E"  w  39  eV. 

At  the  lower  energy  E'  ■  13  eV,  the  calculation 
does  not  match  the  position  of  the  experimental  inter- 
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mediate  peak.  In  this  case  the  diabatic  probability  is 
significantly  leu  than  O.S,  and  (1  -  P)P,  the  overall 
branching  probability  for  following  the  ionic  path, 
decreases  strongly  u  b'  increases  toward  Re.  This 
decreases  the  low -angle  side  of  the  intermediate  peak 
and  causes  the  experimental  peak  position  to  shift  to 
larger  r  values. 

Better  agreement  with  the  experimental  data  at  13 
eV  could  obviously  be  achieved  by  using  estimates  of 
the  coupling  matrix  elements  to  weight  the  differen¬ 
tial  cros  sections  by  Landau-Zener  curve-crossing 
probabilities  f*(l  -  P)  as  a  function  of  b.  An  even 
better  fit  would  be  achieved  by  using  improved  poten¬ 
tials  in  the  crouing  region  [6] .  However,  the  point  of 
the  present  example  is  that  the  existence  of  the  inter¬ 
mediate  peak  in  the  angular  distribution  can  be 
explained  without  invoking  these  complications.  Varia¬ 
tion  of  P  with  b  and  adiabatic  corrections  to  the  poten¬ 
tial  can  be  viewed  as  secondary  effects  that  do  not 
cause  the  intermediate  peak,  but  merely  modify  its 
shape  and  position. 

4.  Discussion 

The  example  in  the  preceding  section  clearly  is  a 
good  system  for  applying  our  simple  model.  The  fact 
that  the  surface  crouing  distance  Rc  is  much  larger 
than  the  region  where  repulsive  terms  become  impor¬ 
tant  implies  a  reasonable  separation  between  the 
secondary  peak  and  the  ionic  rainbow.  This  separation 
is  evident  from  a  comparison  of  the  impact  param¬ 
eters  corresponding  to  the  secondary  peak  [ij  *  6.3  A 
from  eq.  (19)1  with  the  impact  parameters  important 
in  the  rainbow  region  (see  fig.  3).  The  large  value  of 
Rc  also  implies  that  our  idealized  potentials  are 
reasonable  approximations  to  the  true  potentials  for 
distances  important  in  the  secondary  peak  region. 
Finally,  reaction  (21 )  does  not  have  any  of  the  com¬ 
plications  associated  with  molecular  targets. 

In  other  systems  where  the  crossing  distance  Rc  is 
closer  to  the  minimum  in  the  ion-pair  potential,  the 
secondary  peak  and  the  ionic  rainbow  peak  will  be 
closer  together.  In  some  cases  they  may  coalesce  into 
one  single  unresolved  peak.  Likewise,  the  covalent 
peak  may  also  be  attenuated  by  repulsive  interactions 
that  spread  the  scattering  over  a  wider  angular  range. 

Molecular  targets  [  1 1  indeed  have  much  lower 
vertical  electron  affinities  than  I  atoms  and  the  sur¬ 
face  crouings  are  located  at  much  smaller  distances. 


In  addition,  bond  stretching  of  the  molecular  ion  on 
the  ion-pair  surface  can  dramatically  modify  branch¬ 
ing  ratios  [7],  Orientation  effects  are  also  important, 
since  coupling  matrix  elements  can  have  strong  depen¬ 
dences  on  collision  geometry  (8).  Finally,  adiabatic 
prestretching  can  yield  ion-pair  products  for  impact 
parameters  larger  than  the  nominal  Rc  value  [9]. 
These  molecular  effects  taken  together  can  shift, 
modify,  enhance,  or  mask  the  intermediate  peak. 
Other  more  detailed,  theoretical  treatments  [7,9,10], 
relying  heavily  on  numerical  calculation,  have  been 
used  to  investigate  these  molecular  effects.  We  have 
no  reason  to  add  these  complications  to  our  very 
simple  model. 

Our  primary  goal  in  this  paper  has  been  the  devel¬ 
opment  of  a  simplified,  analytic  model  that  illumi¬ 
nates  the  basic  cause  for  the  secondary  peak  some¬ 
times  observed  in  ion-pair  formation  differential 
cross  sections.  Knowledge  of  the  simplest  framework 
for  describing  this  peak  should  aid  in  the  application 
and  interpretation  of  more  sophisticated  calculations. 
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For  Ar*  *  Oj  collisions  at  a  100  eV  energies  sn  intense  rainbow  peak  has  been  observed 
In  the  neutral-product  angular  distribution.  The  surprisingly  large  inelasticities  (1-2  eV) 
and  the  absence  of  elastic  scattering  in  this  small-angle  feature  are  explained  by  a  charge- 
transfer  model  involving  transient  formation  of  an  Ar*  *  Oz '  Intermediate.  Large  pro¬ 
duct  inelasticities  are  understood  by  examing  the  relative  timii^  of  the  projectile  motion 
and  the  Oj '  vibration  on  the  intermediate  surface. 


PACS  numbers:  34.SO.Ea,  34.10,+x,  34. SO 

We  describe  detailed  double-differential  cross- 
section  measurements  for  the  interaction  of  a 
metastable-argon  beam  with  02  target  molecules. 
Ion-pair  surfaces  (Ar*  +02")  play  a  key  role  in 
this  system  and  indeed  we  find  that  Ar*  ions  are 
an  important  reaction  product.  Here  we  present 
dramatic  evidence  that  the  ion-pair  surface  also 
acts  as  an  intermediate  state  in  an  important  in¬ 
elastic  energy  transfer  process  yielding  neutral 
products.  This  process  yields  unexpectedly  large 
values  of  translational-  to  internal -energy  trans¬ 
fer  for  interactions  occurring  at  relatively  large 
Impact  parameters.  The  large  inelasticities  and 
the  apparent  absence  of  elastic  scattering  for 
these  ion-pair-mediated  collisions  will  be  ex¬ 
plained  here  using  a  very  simple  model  that  com¬ 
pares  the  timing  of  the  relative  translational  mo¬ 
tion  to  that  of  the  vibrational  motion  of  the  transi¬ 
ent  02*  intermediate  ion.  In  this  work,  we  are 
probing  the  crucial  energy  range  where  the  colli - 
sional  time  and  vibrational  period  are  of  com¬ 
parable  size. 

Our  apparatus1  and  techniques2  have  been  de¬ 
scribed  previously.  Briefly,  an  Ar beam, 
produced  by  near-resonant  charge  transfer  of 
Ar*  in  Rb  vapor,  enters  a  collision  cell  filled 
with  02.  The  angular-energy  distribution  of  Ar* 
product  ions  from 

Ar*+02-Ar*+02*  (1) 

is  measured  by  a  rotatable  Channeltron  mounted 
behind  a  127°  electrostatic  energy  analyzer.  A 
second  rotatable  Channeltron  detector  is  used  to 
determine  the  angular  distribution  of  scattered 
neutral  products.  Time-of-flight  (TOF)  tech¬ 
niques  are  used  to  measure  velocity  distributions 


S,  34.70.+e 

and  determine  inelasticities  for  these  scattered 
neutrals. 

The  angular  distribution  of  Ar*  product  ions 
seen  in  Fig.  1(a)  has  a  single  peak  at  1.4°,  an  in¬ 
tense  rainbow  feature  caused  by  the  potential 
well  of  the  Ar*  +02"  surface  [see  Fig.  2(a)]  for 
trajectories  that  transfer  to  the  ion-pair  surface 
on  the  initial  inward  passage  through  the  surface 
crossing. 

We  observe  a  similar  rainbow  peak  at  a  nearly 
identical  angle  in  measurements  of  the  neutral- 
product  angular  distribution  [Fig.  1(b)].  This 
rainbow  must  also  result  from  trajectories  that 
sampled  the  ionic  potential  well  during  the  colli¬ 
sion,  since  purely  repulsive  scattering  on  the  up¬ 
per  covalent  potential  would  yield  a  monotonic, 
nearly  flat,  p  vs  6  distribution.  A  second,  near¬ 
ly  flat,  underlying  contribution  evident  in  Fig. 

1(b)  can  be  associated  with  these  purely  repulsive 
interactions. 

TOF  measurements  give  additional  evidence  for 
the  two  separate  major  pathways  to  neutral  prod¬ 
ucts.  Figure  3  shows  a  typical  Ar*  TOF  distribu¬ 
tion  at  an  angle  in  the  rainbow  region.  The  dis¬ 
tribution  can  be  characterized  by  its  nominal 
c.m.  exothermicity  value  Q  (--0.73  eV)  and  its 
energy  width  W  [-  3.0  eV  full  width  at  half  maxi¬ 
mum  (FWHM)].  When  measured  values  of  W  and 
Q  are  plotted  vs  9  in  Figs.  1(c)  and  1(d),  it  is  evi¬ 
dent  that  the  rainbow  region  has  anomalously 
large  values  of  inelasticity  and  inelastic  widths 
compared  to  the  distributions  at  aiuer  angles. 

The  underlying  covalent  repulsive  scattering  has 
a  gradually  increasing  inelasticity  (and  a  slightly 
increasing  width)  as  the  angle  increases.  This  is 
consistent  with  a  simple  momentum-transfer 
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FIG.  1.  Date  for  Ar*  +  Oj  collisions  at  a  beam  en¬ 
ergy  E,  •  114  eV.  (a)  Product  Ar*  laboratory  angular 
distribution;  plotted  la  p  =  8sin Bdo/du.  (b)  Scattered 
Ar*  angular  distribution,  (c)  Half-width  of  the  Ar* 
energy  distribution  vs  0.  (d)  Ar*  c.m.  energy  loss  vs 
laboratory  0. 


> 

V 
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FIG.  2.  Cuts  through  the  two  (idealized)  potential 
surfaces,  (a)  A  cut  at  the  O,  equilibrium  separation 
r=  r,  ,  showing  the  initial  surface  crossing  at  R  »  Rc 
»  3.44  A  .  R  is  the  distance  from  the  atom  to  the  cen- 


process  at  small  H.  The  width,  average  Q  value, 
and  fractional  contribution  of  the  covalent  compo¬ 
nent  can  readily  be  extrapolated  into  the  rainbow 
region,  as  shown  by  the  dashed  lines  in  Fig.  1. 
This  allows  one  to  estimate  and  remove  the  co¬ 
valent  contribution  from  TOF  distributions  meas¬ 
ured  in  the  rainbow  region  and  to  obtain  the  re¬ 
sultant,  highly  inelastic,  contribution  from  the 
ion-pair-mediated  process  (see,  e.g..  Fig.  3), 

Ar*  +02— [Ar**Oa*]  — Ar*  +Ot*Q.  (2) 

Similar  large  inelasticities  are  obtained  at  other 
angles  in  the  rainbow  region  for  process  (2).  At 
lower  energies  (£„-  70-90  eV),  this  inelasticity 
decreases  slightly;  otherwise  the  results  are 
quite  similar  to  those  shown  here.  Clearly  proc¬ 
ess  (2)  yields  very  little  elastic  (Q  *0)  scatter¬ 
ing.  The  true  elastic  contribution  is  even  small¬ 
er  than  it  appears  in  Fig.  3,  since  the  TOF  dis¬ 
tributions  have  not  been  corrected  for  experimen¬ 
tal  resolution  (an  apparatus  resolution  FWHM, 
estimated  from  TOF  beam  profiles,  is  shown  in 
Fig.  3). 


ter  of  the  molecule.  (b>  Cuts  at  fixed  R.  At  the  Initial 
idealized  crossing  (R  =  R^  =  3.44  A)  the  O-,'  is  formed 
compressed;  as  it  vibrates  (changes  of  r  along  hori¬ 
zontal  line  in  O,"  potential  well),  the  surface  crossing 
distance  at  the  second  crossing  Rc’  also  varies  (3.44  A 
*  Re'  *  5.78  A).  The  Ar*  +  Oj  curve  is  essentially  un¬ 
changed  for  3.44  A  R  a  5.78  A ,  while  the  Ar*  +  Oj  * 
curve  moves  vertically  *  -  1/ft . 


We  believe  that  the  large  inelasticities  and  neg¬ 
ligible  elastic  scattering  for  process  (2)  in  the 
rainbow  region  can  be  understood  as  a  simple 
consequence  of  the  timing  of  02*  bond  stretching 
(and  subsequent  vibration)  between  the  two  sur¬ 
face  crossings  on  the  incoming  and  outgoing  por¬ 
tion  of  the  trajectory.  Vibrational  timing1  has  in¬ 
deed  been  invoked4  to  explain  oscillations  in  the 
energy  dependence  of  total  cross  sections  for  ion- 
pair  formation  in  the  similar  Cs  and  K  +0,  sys¬ 
tems.  The  effects  of  O,'  vibrational  timing  have 
also  been  observed  quite  recently  in  the  product 
angular1  and  energy*  distributions  for  these  re¬ 
lated  systems. 
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FIG.  3.  Typical  Ar*  TOF  (and  energy-loss)  spectrum 
in  the  rainbow  region.  The  observed  spectrum  (heavy 
line)  has  contributions  from  (a)  pure  covalent  scattering 
(dashed  curve)  and  (b)  lon-pair— mediated  inelastic 
scattering  (light  solid  line).  The  dashed  curve  is  es¬ 
timated  from  TOF  distributions  obtains  Just  beyond 
the  rainbow  region,  shifted  to  die  energy  loss  at  1.62* 
indicated  by  die  dashed  portion  of  Fig.  1(d),  and  nor¬ 
malized  to  give  branching  ratios  consistent  with  the 
partitioning  in  Fig.  1(b).  The  light  solid  cptvs  is  ob¬ 
tained  by  subtraction. 


We  use  a  simplified  classical  picture  to  illumi¬ 
nate  the  basic  physics  underlying  the  mechanism. 
Assume  that  the  O,'  and  02  potential  curves  are 
not  perturbed  by  Ar*  and  Ar*,  respectively,  at 
distances  R>Re-  3.44  A,  that  the  Coulombic  term 
-  1/R  dominates  the  ion-pair  interaction  at  these 
distances,  that  the  02  begins  at  rest  at  the  bottom 
of  its  well  at  r,(02),  that  transitions  between  dia- 
batic  electronic  states  occur  vertically  (with  r 
unchanged)  and  only  at  the  surface  crossing  seam, 
that  electron  transfer  (at  t  =  0)  initiates  a  vibra¬ 
tion  (initially  an  expansion)  of  02* ,  and  that  the 
motions  along  r  and  R  on  the  ion-pair  surface 
evolve  independently.  Figure  2(b)  shows  the  Ar* 
+02  and  Ar*  +Oa"  surfaces  at  R  =  3.44  A,  the  dis¬ 
tance  (R  *RC)  corresponding  to  the  initial  transi¬ 
tion  to  the  ionic  surface,  in  our  simple  picture, 
the  02*  is  produced  with  a  particular  internal  en¬ 
ergy  (shown  by  the  horizontal  line  going  to  the 
right  from  the  intersection  point).  As  R  decreas¬ 
es  to  its  minimum  value  (at  the  turning  point)  and 
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FIG.  4.  Dependence  of  Rc'  (solid  lines)  and  R  (dashed 
lines)  on  tune  t  ?fter  the  initial  transfer  to  the  ion- 
pair  surface.  Rc '  is  shown  (unlabeled)  for  the  incoming 
neutral  electronic  states  and  (labeled)  for  several 
higher  neutral  states  of  Ar  or  Oj .  RV  ,b,E0)  is  shown 
for  several  6  (in  A)  and  for  E„  »  114  eV. 


then  increases,  the  02*  vibrates.  A  knowledge  of 
the  shapes  of  the  02  and  Oz*  curves  allows  an  es¬ 
timate4  of  the  time  dependence  of  the  surface 
crossing  radius  Rc'(t)  to  be  made  from  the  classi¬ 
cal  vibrational  timing,  r(f ),  of  the  O/.  Rc'ft)  de¬ 
pends  on  the  Ar*  ionization  potential,  <p,  f ,  and 
the  time -dependent  vertical  electron  affinity  of 
02,  (PeaM*)),  via  Re’(f)<=  (14.4  eV  A)/[<Pn>- <Pea 
x  (rft ))].  The  relative  timing  of  the  translational 
motion  and  the  02*  vibration  determine  the  r  and 
R  values  for  the  second  surface  crossing.  With 
proper  timing  a  fully  expanded  02‘  (after  0.5  vi¬ 
brations)  could  recross  to  the  covale.it  surface  at 
R  =5.78  A  [see  Fig.  2(b)],  producing  an  inelastic 
excitation  of  Oa  by  ~  1.7  eV. 

Figure  4  shows  the  dependence  of  Re'  and  R  on 
the  time  t  after  the  initial  surface  crossing  on 
the  incoming  trajectory.  Rit)  was  calculated  us¬ 
ing  straight-line  trajectories  for  simplicity.  For 
a  given  b,  the  intersection  of  R(t)  and  R/(t)  de¬ 
termines  the  time  of  the  second  crossing  and  the 
inelasticity  of  reaction  (2),  Q  =  - 14.4 (1/RC  - 1/ 
Rt'),  for  those  collisions  returning  to  the  cova¬ 
lent  surface  at  the  second  crossing.  Note  that 
the  initial  expansion  of  Re'  is  so  rapid  that  tra¬ 
jectories  with  b  values  Just  inside  of  the  cross¬ 
ing  at  R  *  3.44  A  cannot  overtake  it  and  must  meet 
it  later  as  it  contracts.  The  long  time  spent  on 
the  ionic  surface  and  the  resultant  significant  at¬ 
tractive  interactions  are  responsible  for  the  rain¬ 
bow  feature  in  the  neutral  scattering  [Fig.  1(b)] 
occurring  at  an  angle  comparable  to  the  rainbow 
peak  angle  in  the  ionic  scattering  [Fig.  1(a)]  of  re- 
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action  (1).  Classical  deflection  function  calcula¬ 
tions7  for  the  ion-pair-mediated  process  (2)  show 
that  the  impact -parameter  range  2.5  A  s  6  s  3.44 
A  contributes  to  the  rainbow  feature.  Hence,  in 
this  simple  picture,  the  rainbow  feature  of  the 
neutrals  would  yield  inelasticities  between  -  1.0 
and  1.6  eV. 

This  simple  model  explains  the  lack  of  elastic 
scattering  in  the  rainbow  region  for  reaction  (2); 
but  it  slightly  underestimates  the  average  inelas¬ 
tic  scattering  values.  Although  relaxing  some  of 
our  simplifying  assumptions  might  yield  a  better 
fit  to  the  data,  the  basic  physics  is  clearly  al¬ 
ready  contained  in  the  model.  One  possible  route 
to  larger  inelasticities  involves  transfer  from  the 
ion-pair  intermediate  surface  to  more  highly  ex¬ 
cited  neutral  surfaces.  The  crossing  radii  Re' 
between  the  ion-pair  surface  and  higher  neutral 
surfaces  also  oscillate  with  time.  A  few  exam¬ 
ples  are  shown  in  the  upper  portion  of  Fig.  4.  Al¬ 
though  the  larger  values  of  R,'  seem  to  imply 
smaller  coupling  matrix  elements  and  thereby 
much  lower  adiabatic  crossing  probabilities,3 
there  may  be  other  factors  favoring  production 
of  electronically  excited  states.  In  fact,  the  radi¬ 
al  motion  of  the  particles  [R «)]  relative  to  that  of 
the  moving  crossing  [rc'(0]  should  be  the  radial 
velocity  appropriate  for  estimating  curve-cross¬ 


ing  probabilities.  In  the  extreme  case,  when  the 
crossing  sphere  is  expanding  rapidly  and  parallel 
to  R(f),  e.g.,  for  6  =  2.5  and  the  2 p,  state  of  Ar  in 
Fig.  4,  the  interaction  can  be  profoundly  altered. 
Significant  electronic  excitation  might  result 
This  work  received  financial  support  from  the 
U.  S.  Office  of  Naval  Research.  Valuable  discus¬ 
sions  with  A.  P.  Hickman,  J.  Los,  and  A.  W. 
Kleyn  are  gratefully  acknowledged. 
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SUBPICOSECOND  EXPERIMENTS  ON  MOLECULAR  VIBRATION 
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Vibrational  motion  on  the  subpicosecond  time  scale  of  a  transient  collision 
intermediate  has  been  observed  through  its  striking  effects  upon  differential  and 
total  cross  sections  in  inelastic  and  reactive  processes  in  atom-  and  ion- 
molecule  collisions.  Transient  charge  transfer  processes  that  form  a  short-lived 
diatomic  Ion  and  then  return  to  the  original  potential  energy  surface  are  shown 
to  be  efficient  paths  to  high  vibrational  excitation  of  molecular  targets* 
Similar  effects  should  be  observable  In  a  large  number  of  collision  systems* 

Introduction 


A  precise  knowledge  of  the  potential  energy  surface  for  a  trlatomlc  system  could. 
In  principle,  be  used  to  predict  the  detailed  distribution  of  product  states  for 
inelastic  or  reactive  processes,  that  probe  the  trlatomlc  surface  as  a  transient 
intermediate  during  a  collision.  in  contrast,  even  when  detailed  atate-to-state 
rates  are  available,  it  is  far  more  difficult  to  extract  quantitative  information 
on  the  relevant  potential  surfaces  from  measurements  of  inelastic  and  reactive 
processes.  The  critical  significance  of  the  collision  Intermediate  Justifies 
strenuous  experimental  efforts  to  improve  our  view  of  its  fleeting  existence* 

One  approach  to  examining  a  collision  intermediate  is  to  prepare  it  directly  in  a 
well-defined  state  or  energy  level  and  observe  its  subsequent  evolution* 
Examples  using  photoexcltatlon  of  neutrals  and  ions  are  given  in  the  chapters  by 
A*  Gluati-Suzor'*'  and  P.  C.  Cosby.  ^  These  "half-collision"  experiments  can 
often  analyse  decay  channels  and  surface  couplings  with  great  precision  assuming 
sufficient  spectroscopic  information  is  available. 

An  alternate  and  quite  difficult  approach  involves  the  examination  of  a  transient 
collision  intermediate  by  direct  spectroscopic  observation  of  the  very  short¬ 
lived  colliding  system  (see,  e.g. ,  chapter  by  J*  Polanyl^)  or  by  optical  pump¬ 
ing  of  the  system  while  it  Is  reacting  (see,  e.g.,  chapter  by  A.  Gallagher^). 
These  experiments,  closely  related  to  spectroscopic  line  broadening,  are  still  in 
an  exploratory  stage  and  must  eventually  address  the  difficult  problem  of  deter¬ 
mining  the  dependence  of  emission  line  widths  upon  nuclear  coordinates  in  the 
transient  intermediate* 

In  contrast,  our  approach  to  examining  the  collision  intermediate  uses  more  tra¬ 
ditional  scattering  techniques,  but  adjusts  experimental  conditions  so  that  the 
transient  Intermediate  Itself  provides  an  internal  clock  wfc-jam  timing  can  be 
observed  through  its  profound  effects  upon  the  resultant  product  distributions* 

^Visiting  scientist  1980-81  at  IBM  Research,  San  Jose,  CA  95193,  U.S.A. 

* Research  support  from  NSF  and  ONR  are  gratefully  acknowledged. 
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The  Internal  clock  Is  the  vibrational  motion  of  a  transient  Ion,  formed  as  part 
of  the  short-lived  intermediate*  We  %rlll  demonstrate  for  a  few  relatively  simple 
systems  Involving  charge  transfer  f  iteraedlate  states  that  very  detailed  informa¬ 
tion  about  the  collision  dynamics  can  be  obtained.  In  some  cases  vibrations  of 
the  transient  Intermediate  on  the  subpicosecond  time  scale  can  be  seen  experi¬ 
mentally.  The  vibrational  notion  of  this  transient  species  yields  the  most 
dramatic  results  in  carefully  chosen  systems  when  the  collision  velocity  is 
adjusted  to  make  the  collision  time  comparable  to  the  vibrational  time. 

Atom-Atom  Ion-Pair  Formation 


Most  of  the  work  described  in  this  report  will  Involve  neutral  atom-molecu?  e  col¬ 
lisions  that  are  strongly  influenced  by  charge  transfer  to  ion-pair  states  during 
the  collision  process*  It  is  therefore  useful  first  to  consider  s  related,  but 
simpler,  atom-atom  Interaction,  that  of  Na  with  1*  Shown  in  Figure  1  are  the 
lowest  l  diabetic  covalent  and  ionic  potential  curves  for  Mai.  At  large  dis¬ 
tances  A,  the  lou-palr  state  is  asymptotically  higher  than  the  covalent  state  by 
&E  *  IP(Na)  -  £A(I),  where  IP  Is  the  ionisation  potential  of  Na  and  £A  la  the 
electron  affinity  of  the  I  atom.  The  ion-pair  state  decreeses  in  energy  as  R 
decreases  and  crosses  the  covalent  state  at  R  ■  RQ.  In  the  simplest  approxima¬ 
tion,  considering  only  the  coulomblc  term  of  V^on,  R£  -  14.4/aE,  with  AE  In  eV 
and  Rc  in  A. 
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Figure  1.  Diabatic  Na-I  potentials. 


Since  the  l£+  covslent  curve  Is  the  same  symmetry  as  the  lon-palr  curve,  there 
will*  be  an  avoidsd  crossing  In  the  corresponding  adiabatic  curvaa  ft+R  -  Rc*  For 
thermal  energy  Na  and  I  atoms  approaching  aach  othar  along  tha  I  curve  (1/8 
of  the  collisions  statistically),  an  aiectron  will  juap  from  Na  to  I  in  thm 
vicinity  of  Rc  as  tha  adiabatic  atats  abruptly  changes  from  a  covalent  to  an 
ionic  charactar.  After  the  Inner  Ionic  region  Is  visited,  tha  outgoing  trajec¬ 
tory  will  once  again  cross  Rc  and  tha  alactron  will  raturn  to  the  Na  . 
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Am  one  raises  che  collision  energy,  one  expects  eventually  that  a  fraction  of  the 
collisions  will  proceed  dlabaclcally  through  the  crossing  region  at  R^  and  that 
for  collision  energies  Eq  >  AE,  the  lon-pair  product  channel  N*+  +  I  would  be 
observed.  Formation  of  ions  indeed  has  been  observed  In  both  total  and  differ¬ 
ential  cross  section  experiments. 


The  probability  for  ion-pair  format ion  depends,  of  course,  on  the  strength  of  the 
coupling  between  the  two  interacting  surfaces.  In  the  simplest  Landau- Zener 
(L-Z)  description  (see,  for  example,  refs.  7,  8),  the  probability  P  for  dlabatlc 
behavior  at  a  crossing  through  the  sphere  of  radius  Rc  is 

2xH? 

1  ‘  (1) 


where  Is  the  coupling  matrix  element,  vR  is  the  radial  velocity  at  R  -  Rc, 

and  $  •"_ / v  -  V  )  evaluated  at  R  ■  R_.  Feist  and  Levine^  have  obtained 

dRv  cov  ion  c 

excellent  fits  to  che  Na  +  I  scattering  data  by  using  Improved  adiabatic  poten- 
tlals  and  a  slight  Modification  of  tha  crossing  distance.  Nonetheless  the  staple 
dlabatlc  IrZ  tsodel  la  a  reasonable  first  approximation  to  the  scattering  process. 


Atom-Molecule  Systems;  Oi  Target 


Collisions  of  fast  alkali  atoaa  with  02  molecules  are  Influenced  In  an  analogous 
way  to  Figure  1  by  the  crossing  between  the  covalent  Incoming  potential  and 
strongly  attractive  Ionic  potentials.  Necessarily  the  three  atoa  system  Is 
vastly  more  complicated  by  the  extra  degrees  of  freedom,  and  the  greater  uncer¬ 
tainties  In  the  knowledge  of  tha  potential  surfaces  and  their  Interactions.  The 
Simple  potential  curves  of  Na  +  t  are  replaced  by  multidimensional  potential 
surfaces  In  the  trlatomlc  system.  Nonetheless,  a  carefully  chosen  set  of 
physically  reasonable  simplifications  allows  us  to  generate  an  Interesting 
classical  picture  of  the  collision  dynamics. 


The  basic  assumption  Is  an  uncoupling  of  the  Internal  notion  of  Che  target  molec¬ 
ule  from  the  relative  translational  motion.  To  apply  this  assumption,  •*  must 
exclude  Impact  parameters  that  explore  the  Inner  repulsive  core  of  the  Inter¬ 
action  potential.  This  neatly  excludes  the  collisions  resulting  In  the  normal 
Impulsive  mechanisms  for  Inelastic  excitation  transfer. 

Consider  a  Cs  atom  approaching  an  02  molecule  on  a  covalent  surface.  Assume  that 
che  0-0  distance  r  Is  the  classical  equlllbrlixi  separation  of  02,  r%.  At  a  Cs-02 
distance,  Rc,  corresponding  to  a  crossing  with  the  Ionic  surface,  an  electron  can 
jump  to  form  Cs*  +  02~.  In  the  spirit  of  our  simple  picture,  sssume  s  vertical 
transition,  an  Instantaneous  electron  Jump  that  forms  02~  with  an  unchanged  bond 
lengch,  r  .  Clearly  then,  the  crossing  distance  Rc  Is  given  by  R£  ■  14.A/[IP(Cs) 
-  EAyCOj)]  where  EAy  Is  the  vertical  electron  affinity  of  02  at  r  -  r#.  Since 
the  equilibria  lnternuclear  separation  of  02~  Is  larger  than  that  of  02,  che 
nascent  02~  begins  its  existence  In  a  compressed  state  on  Its  repulsive  lnnar 
wall  (see  che  02,  02~  potentials  in  Figure  2).  As  tha  Cs*  and  02~  Iona  axplore 
che  attractive  region  Inside  Rc,  the  02~  starts  to  expand,  thus  Initiating  an  02~ 
vibration.  It  Is  clear  from  Figure  2  that  as  Che  02“  goes  through  one  complete 
classical  vibration,  the  vertical  electron  affinity  (which  depends  on  the  Instan¬ 
taneous  value  of  r)  first  Increases,  then  decreases  with  a  period  Identical  to 
that  of  the  02”  vibration.  Since  Re  depends  on  the  Instantaneous  value  of  EAy(r) 
through 


Rc  -  l*.A/(IF(Ca)  -  EAvlr0j-(t)l),  (2) 
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this  crossing  radius  oust  vibrate  In  and  out  with  the  seas  period  aa  the  Oj*. 
Hence  the  position  of  the  second  crossing  with  the  covalent  surfece  depends 
critically  on  the  relative  tlalng  of  the  translational  notion  and  the  Oj'  vibra¬ 
tion.  Since  che  coupling  decreases  exponentially  with  the  branch¬ 
ing  probability  P  at  thla  second  crossing  Is  also  a  very  aenaltlve  function  of 
che  relative  tlalng.  Depending  on  collision  velocity  and  lapact  paraaeter,  the 
transient  lon-palr  lntemediate  aay  have  a  second  crossing  with  the  covalent 
surface  either  In  an  expanded  geoaetry  that  favors  diabetic  behavior  or  In  a 
coapreaaed  geoaetry  chat  favors  adiabatic  behavior  In  che  coupling  region. 


Figure  3  shows  scheaatlcally  Cs  +  Oj  collisions  ac  -  35. S  eV  and  various 
lapact  paraaetera.  Each  trajectory  (drawn  as  a  straight  line  for  alapllclty  of 
presentation)  la  assumed  to  convert  to  Ce+  +  0j“  lon-palre  et  the  Initial 
crossing  at  R  -  R(.  At  that  Instant  (t  *  0)  the  crossing  radius  »c  starts 
vibrating  with  the  Oj“  period  In  accord  with  Equation  (2).  The  calculated  radius 
Eg  et  the  second  crossing  point  (a  function  Of  b)  la  Indicated  In  Figure  3.  The 
lower  trajectories  labelled  l,  2,  3  have  lapact  parses  ten  that  correspond  to 
ranching  the  second  crossing  after  1,  2,  or  3  full  vibrations  of  che  02~  (and  of 
Eg).  These  trajectories  would  tend  to  favor  adiabatic  passage  through  the  second 
crossing  region  yielding  raneutrallxed  products.  In  contreet  the  upper  trajec¬ 
tories  have  lapact  paraaacers  that  would  reach  the  second  crossing  after  -  1.5 
or  2.5  Oj"  vlbratlona  and  would  yield  larger  diabetic  probabilities  for  lon-palr 
products.  Clearly  the  branching  ratio  would  have  en  lapact  paraaeter  dependence 
that  alght  show  up  as  oscillations  In  the  differential  cross  sections  for  this 
syscaa. 


Ci»Oj  JS.SsV 


Figure  3.  Vibrational  tlalng 
effecca  on  R£  for  Cs  +  Oj 
collisions.  The  straight  line 
trajectories  (et  boctoa)  favor¬ 
ing  re neutral Isa cion  would 
suffer  quite  different  attrac¬ 
tive  deflections  during  che  else 
spent  on  the  ionic  lntemediate 
surface.  The  Insert  at  the 
right  indicates  angular  posi¬ 
tions  end  relative  intensities 
for  three  reneutrallaatloe  peeks 
predicted  frea  classical  trajec¬ 
tory  calculations. 
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Figure  4  shows  chs  angular  distribution  ot  osutrsl  scattering  in  vsrious  colli¬ 
sion  systsms.111)  The  35.3  eV  Cs  +  Oj  date  Indeed  has  distinct  structure  that  is 
conalstent  with  the  suggested  vibrational  timing  arguaents  in  Figure  3.  This  can 
be  coapared  to  the  nearly  flat  distribution,  associated  with  purely  repulsive 
potentials,  evident  for  Cs  *  Ar  scattering  also  in  Fig.  4.  At  higher  collision 
velocities,  the  structure  for  Cs  +  Oj  changes  draastically  since  the  reneutral- 
lxlng  trajectories  shift  to  smaller  impact  parameters  and  the  number  of  observ¬ 
able  vibrational  oscillations  decreases  as  the  collision  time  decreases.  The 
value  of  R£  not  only  influences  the  branching  probability  for  reneutralixatlon, 
but  also  strongly  affects  the  extent  of  deflection  since  R£  is  the  distance  at 
which  the  coulombic  force  is  'switched  off*;  its  variation  implies  a  related 
veriation  in  the  total  time  that  the  attractive  force  operates. 


Figure  4.  Neutral  product  angular 
distributions  plotted  vs.  reduces 
angle  x  •  M  for  several  collision 
systsas  at  various  beam  energies. 


The  energy  dependence  of  the  total 
cross  section  for  ion-pair  formation 
in  the  Cs  02  system1 123  is  shown 
in  Fig.  5.  Although  integrating  the 
oscillating  ion-pair  fotmatlon 
probability  over  all  lnpaet 
parameters  at  each  energy  might  be 
expected  to  wash  out  any  strong 
oscillations  in  the  velocity 
dependence  of  the  total  cross 
section,  one  still  observes  these 
oscillations.  In  fact,  the  inter¬ 
pretation  of  oscillatory  total  cross 
section  data  for  ion-pair  formation 
in  Cs  +  Oj  and  K  +  Oj  collisions113) 
yielded  the  first  dramatic 
indication  that  miltlple  transient 
molecular  vibrations  could  be  used 
J  as  an  Internal  dock  for  the 
collision  intermediate. 

Returning  to  the  differential  cross 
'  section  measurements  In  Figure  4,  we 
note  that  the  lighter  K  atom 

explores  a  higher  velocity  regime 

than  Cs  in  the  same  energy  range. 
In  fact,  the  reneutralixatlon 

structure  simplifies  considerably 
for  K  *  Oj  at  90  eV.  In  this  case 
the  time  spent  in  the  inner  ionic 
region  is  somewhat  shorter  than  the 
time  for  one  0j~  vibration  for  all 
impact  parameters  larger  than  the 
extent  of  the  repulsive  inner 

core.  The  single  peak  is  a  rainbow 
feature  in  the  reneutrallsatlon 
angular  distribution.  Figure  i 
shows  a  similar  reneutrallsatlon 
rainbow  in  the  scattering  of 
matesCable  Ar  atoms  by  Oj  at 
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Figure  5.  Total  Ionization  cross 
section  for  C*  +  Oj  collisions. 
Closed  circles  show  negative  Ion 
products;  open  circles— electrons. 


Figure  6.  Ar  +  0,  ecatterlng 
*t  !,  +  114  eV.  (a)  Product 
Ar  angular  distribution 
p  -  S^sin  8  do/ du  vs.  8. 

(b)  Ar  neutral  angular 
distribution,  (c}  Half -width 
of  the  product  Ar  energy 
profile  vs.  8.  (d)  Average 

Ar  c.  is.  energy  loss  vs. 
laboratory  8. 


E  •  114  •».<»>,  at  +  0j+-»  Ar  +02  ♦  Ar  0j.  The  corresponding  lon-palr 
rainbow  (or  Ar  +  02  ♦  Ar  +  0j  Is  also  shown.  The  neutral  rainbow  Is  located 
at  slightly  saallsr  angles  because  the  couloablc  attraction  Is  cut  off  at  IJ  when 
rcneutrallsatlon  occurs.  The  slallarltles  of  the  Ar*  and  K  reactions  with  (>2  are 
not  surprising  and  ere  consistent  with  an  nr  other  coaparlson  studies  between 
ae castable  rare  gas  and  alkali  reactions. 


Another  iaportant  consequence  of  the  vibrational  tlalng  upon  the  reneutrallsatlon 
products  la  the  very  large  vibrational  excitation  of  0j  that  can  tesult.  A 
scratched  02*  would  reneutralise  vertically  to  a  vlbratlonally  excited  O2  as  can 
be  seen  in  Figure  2. 
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Indeed,  large  vibrational  Oj  excitation  for  acatterlng  at  anglee  In  the  rainbow 
region  la  observed  b y  clae-of- flight  (TOP)  aeasureaents  of  acattered  Ar*'1  '  (and 
la  the  analogoua  K  +  Oj  eystaa^ *“•***).  When  the  TOP  Inelasticity  Q  and  Its 
■assured  half-width  i  are  plotted  vs.  8  (Pig.  8),  an  Isolated  hlghly-lnelastlc 
feature  appears  prectaaly  In  tha  rainbow  region  built  upon  a  second  feature  whoee 
Inelasticity  slowly  rises  with  8.  The  second  feature  la  associated  with  the 
■ore  coaaon  Impulsive  Inelastic  scattering  on  the  repulsive  inner  wall  of  the 
covalent  potential.  This  covalent  contribution  can  be  seen  in  the  nearly  flat 
background  that  underlies  the  sharp  rainbow  peak  in  the  neutral  angular 
distribution. 

Atoa-Holecula  System  Halogen  Targets 

In  a  saries  of  experlaants  at  tha  POM  Institute,  reported  at  the  Paris  ICPEAC^20^ 
and  reviewed  elsewhere, the  effects  on  lon-pair  fornatlon  of  bond  stretch- 
log  la  alkali  atoa  halogen  nolecule  collisions  were  studied.  The  results  have 
since  been  confirmed  in  work  by  other  groups.  Reactions  of  alkalis  and 
Ar  with  halogens  have  analogous  Mechanises  to  those  described  for  Oj  targets. 
However,  there  are  great  differences  In  Che  details  of  the  potential  aurfaces  and 
tha  resulting  product  distributions.  Plrst,  since  the  vibration  times  of  the 
heavy  halogen  negative  Iona  are  a  factor  of  5-10  larger  than  that  of  Oj",  Che 
halogen  negative  Ion  will  only  undergo  a  fraction  of  a  vibration  befora  the 
second  crossing  le  reached.  Second,  the  larger  vertical  electron  affinity  of  the 
halogens  and  their  mich  more  rapid  Increase  with  r  Imply  a  larger  Rc  and  a  very 
rapidly  Increasing  tc(t).  Hence  In  aany  eases,  an  adiabatic  first  passage  through 
Rc  Will  yield  such  a  rapid  Increase  with  t  of  kc  that  tha  second  crossing  will  be 
alaost  entirely  diabetic.  Hence  tha  ton-pair  rainbow  (as  In  Pig.  6)  will. have  no 
accompanying  neutral  rainbow  froa  reneutralltaclon  of  tha  second  crossing. 

Tet  strong  rainbows  In  the  neutral  scattering  have  been  obeerved  and  have  been 
attributed  to  another  consequence  of  bond  stretching,  vlbratlonally-lnduced  third 
crossings.  Consider  the  covalent  route  to  the  ionic  surface  with  electron  trans¬ 
fer  occurring  at  tha  second  crossing  of  the  sphere  of  radius  *c.  As  the  lone 
continue  to  larger  distances,  the  Brj”  bond  will  start  to  stretch  and  *c  will 
start  to  follow  tha  departing  Ions.  Sines  the  ultimate  Increase  In  k£  for  halo¬ 
gens  Is  very  large,  the  expanding  crossing  radius  will  soon  overtake  tha  outgoing 
K+  Ion.  This  results  In  a  third  crossing  created  by  the  lntraaolecular  notion  of 
the  Brj-.  Tha  situation  Is  shown  schematically  in  Pig.  7  for  K  +  Br2  at  20  eV. 


Plgure  7.  Simplified  picture  of 
a  collision  of  a  fast  K  aton 
and  a  fixed  Sr2  molecule.  The 
repulsive  part  of  tha  potential 
Is  Indicated  by  tha  sphere. 

Two  types  of  trajectories 
leading  to  neutral  products  are 
shown.  Whenever  the  system  is 
neutral  an  0  la  Inserted  In  the 
trajectory;  tha  Ionised  state 
Is  Indicated  with  a  +. 


K  ♦  8rj  20  tV  _ _ 

o _ 2 _ 2. 
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It  re neutralisation  occur*,  r  >  r^Srj)  and  consequently  vibrational  excitation 
will  taka  place.  Thla  vibrational  excitation  indeed  baa  been  obeerved  by  Kleyn 
et  al.'24^  Since  the  couloab  force  1*  only  active  for  a  very  ahort  tie*  at  the 
end  of  the  colllaion  the  scattering  angles  are  saall.  Thus  the  first  effect  of  a 
third  crosalng  occurring  for  covalent  reneutrallsatlon  la  that  anoaalous  vibra¬ 
tional  excitation  can  be  observed  at  very  snail  scattering  angle*.  The  second 
effect  of  covalent  reneutrallsatlon  Is  that  It  leads  to  a  vlbratlonally  Induced 
rainbow.  The  origin  of  this  rainbow  can  be  seen  by  referring  to  Fig.  7.  Ac 
t  •  0,  Vc  starts  Increasing,  eventually  overtaking  the  departing  ton.  For  large 
b  (near  Rc),  the  Ion  trajectory  has  a  vary  snail  coaponent  In  the  R  direction.  Is 
overtaken  quickly  by  the  expanding  sphere,  and  suffers  little  attractive  deflec¬ 
tion.  As  b  decreases,  the  tine  spent  on  the  coulonblc  surface  (before  the  third 
reneutrallxlng  crossing)  increases.  However,  the  resulting  Increase  In  the 
attractive  interaction  with  decreasing  b  Is  eventually  coapansated  by  the  fact 
that  the  attractive  force  has  an  aver  tnaller  coaponent  perpendicular  to  the 
trajectory.  Hence  a  aaxlaua  attraction  can  occur  in  the  deflection  function  for 
lapact  paraaetars  larger  than  Che  repulsive  core.  This  yields  a  rainbow  that  is 
shaped  aore  by  Che  bond-stretching  dynaalcs  than  by  the  R  dependence  of  the  lon- 
palr  potential.  This  rainbow  feature  occurs  at  very  sasll  scattering  angles.  An 
example  (Fig.  8)  Is  shown  for  collisions  of  Ar*  with  I2  at  60  *V. ^*)  The  cross 
section  at  angles  beyond  the  rainbow  region  Is  flat,  and  due  to  underlying 
contributions  froa  repulsive  scattering.  Slailar  behavior  has  been  found  for 
K  +  Brj  collisions.24 


Figure  8.  Neutral  product 
angular  distribution 
p  •  9  sin  8  -7—  plotted  vs. 
reduced  angla““t  -  E.8  for 
Ar  +  at  61.3  *V. 


40  80  120 

r  (eV-degreetl 


Theoretical  Hods  Is  and  Extension  to  Other  Syste 


The  staple  pictorial  aodels  described  here  csn  be  eade  aore  quantitative  by  using 
classical  trajectory  surface  hopping  calculations  on  dlsbatlc  potentials,  idille 
still  assualng  an  uncoupling  between  the  r  and  R  notions  on  Che  two  potential 
surfaces.  The  calculations!  techniques  and  chair  application  to  lon-pelr  forma¬ 
tion  in  collisions  of  alkali  atone  and  Ar  with  'halogen  aolecule*  have  been 
described  in  detail  (see,  e.g. ,  8,  17,  24,  24).  These  calculations!  techniques 
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have  bno  extended  Co  eh*  ijticmi  considered  here,  both  In  their  simplest  font 
and  »lth  extensions  to  simulate  adiabatic  affacte  laportant  for  lapect  parameters 
near  R^.  All  of  the  laportant  observed  features — structure  in  the  neutral 
angular  distributions  with  Oj  targets,  large  Oj  excitations  asdlated  by  the  ion- 
pair  surface,  oscillations  in  the  total  ion-pair  format too  cross  section,  cova¬ 
lent  reneutralisation  rainbows  la  the  scattering  by  halogen  aolecules— appear  in 
the  calculations. 

The  discussion  of  vibrational  notion  so  fsr  has  been  completely  classical.  This 
approach  gives  saclsfectory  qualitative  explanations  for  the  processes  described 
above.  If  one  is  Interested  in  the  description  of  final  vibrational  state  dis¬ 
tributions  obtained  using  stats  specific  dstsctlon,  better  methods  are  required. 
One  method,  the  moving  umvepacket  approach  Introduced  by  Clalason  et  al. 
has  bean  very  successful  in  explaining  final  vibrational  state  distributions  from 
quenching  experiments  and  is  similar  la  spirit  to  the  work  discussed  so  far.  In 
the  moving  wavepackat  method,  the  tine  evolution  of  the  vibrational  wavefuactlon 
Is  represented  by  the  evolution  of  s  wavepackat  on  the  respective  potential  sur¬ 
faces.  When  the  potential  surfacea  are  assumed  additive  in  R  and  r  the  method  is 
easy  to  use.  The  easiest  and  anthanatlcally  most  elegant  version  of  the  eat  hod 
employs  diabetic  potential  surfacea  and  the  wavepackat  is  projected  from  one 
molecular  potential  to  another  at  each  adiabatic  passing  of  a  crossing.  The 
expression  for  the  final  vibrational  state  population  can  be  given  in  cloaed 
fora. 

The  diabetic  modal  has  bean  used  to  explain  energy-dependent  oeclllatlons  in  the 
final  vibrational  •cate  distribution  of  N2(C3rt  )  excited  in  the  quenching  of  Ar* 
by  ground  stare  Hj.  CJ7  Hie  results  are  stioen  in  Fig.  9.  The  experimental 

ratio  (v  -  1/v  -  0)  oaclllates  as  a  function  of  energy.  The  model  calculatlona 
show  that  the  oscillations  sre  due  to  several  vibrations  of  the  intermediate  n2 
molecule,  which  obtained  some  N2~  Character  from  a  higher  Ar+  +  Nj"  surface.  The 
molecular  vibrations  during  the  collisions  srs  reflectsd  in  the  oscillatory 
motion  of  the  wavepackat  in  its  potential.  The  adiabatic  version  of  the  moving 
wavepackat  model  haa  been  used  to  explain  quenching  of  excited  Ha  by  N,  and  CO 
and  very  good  agreement  between  theory  and  experiments  has  been  obtained.^*8' 


Figure  9.  Population  rat  Jo  . 
(v’-O/v’-l)  for  the  H  (C  D  ) 
product  resulting  from*  u 
quenching  of  Ar  by  N,  as  a 
function  of  energy.  The 
circles  sre  the  experimental 
points;  the  solid  line  is  the 
theoretical  result.  The  number 
N  refers  to  the  number  of  half 
vlbraclons  of  the  transient  H2 
molecules. 


The  moving  wavepackat  nodal  is  an  Intermediate  between  the  simple  classics? 
approach  mainly  used  to  explain  this  kind  of  atom-molecule  collisions  and  rare 
complete  quantal  studies.  Klo^  and  Los'1  '  solved  the  eloea  coupled  equations 
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tor  some  Halting  cists  and  daaonatratad  tan  potential  Importance  of  oon-vartlcal 
craoaltlooa  ac  lots  aoarglad.  Rickman/3^  using  a  classical  path  asthod  and  th< 
Magnus  approximation,  observed  oscillations  In  hla  calculated  cross  section  for 
lon-pair  foraatloo  In  Cs  +  02  collisions.  beastly,  bekar  and  Sexoo'313  par- 
forms  d  close  coupling  calculations  for  loalaatic  t  +  Oj  collisions  Involving  the 
lon-pair  Intermediate.  These  calculations  conf In  the  staple  Ideas  developed 
using  the  classical  nodels.  In  addition,  this  kind  of  calculations  gives  such 
more  insight  into  the  detailed  collision  dynamics  and  serves  as  a  standard  with 
which  the  classical  and  seoiclasslcal  techniques  can  be  coopered  and  Judged. 

The  effects  of  vibrational  timing  in  the  subplcosecoud  range  have  been  nalnly 
studied  systaosclcally  In  collisions  of  alkali  or  ontaecable  race  gas  atoms  with 
simple  molecules  lib  halogens,  O2  and  Mj.  There  Is  no  reason  to  believe  that 
these  effects  will  not  occur  In  many  other  collision  systems.  One  obvious  exten¬ 
sion  Is  to  the  study  of  lon-pair  mediated  effects  In  albll  collisions  with 
larger  molecules,  bcently,  the  Importance  of  'bond-bending*  effects  In  the 
Ion'3*'  have  been  discussed  In  analogy  to  bond  stretching  In  the  halogens.  Other 
recent  experiments  have  indicated  the  possible  Involvement  of  similar  vibration- 
alls  inelastic  effects  due  to  temporary  charge  transfer  from  H+  or  H~  to 
02.  33,3*'  Colllslonally-producsd  amissions  In  metastable  rare  gas-halogen  reec- 
clons  at  near  thermal  energies  seem  to  be  similarly  unaffected  by  bond  stretch¬ 
ing.  (35)  In  addition,  these  effects  should  be  obeervable  In  many  excitation 
transfer  or  Penning  ionisation  reactions. 


Ua  have  emphasized  here  the  importance  In  atom-molecule  collisions  of  vibrational 
excitation  that  is  Induced  by  Involvement  of  another  electronic  state  of  the 
transient  Intermediate  species.  In  contrast  to  the  direct  excitation  of  a 
vibrator  by  repulsive  Interactions  on  a  single  potential  surface,  this  Inter¬ 
mediate-Induced  excitation  can  yield  high  Internal  energies  In  large  Impact 
parameter  collisions.  Vibrational  timing  of  the  Intermediate  species  on  a 
subpicosecond  tins  scale  can  be  observed  as  structure  In  neasuraments  of  both 
differential  and  total  cross  sections. 

There  are  several  general  requirements  for  direct  observation  of  these  effects. 
First,  Che  collision  must  Involve  another  electronic  state  of  Che  diatomic  molec¬ 
ule  having  a  somewhat  different  geometry  from  that  of  Che  Initial  state.  In  an 
adiabatic  description,  the  intermediate  must  access  a  region  of  the  potential 
surface,  where  Che  nolecular  bond-strength  Is  drastically  modified  from  chat  of 
the  Isolated  molecule.  System  with  Important  low-lying  lon-pair  states  satisfy 
this  requirement.  The  second  requirement  for  observation  of  vibrational  timing 
la  that  the  transfer  from  one  dlabadc  state  to  the  other  occurs  In  a  region  of 
the  potential  surface  where  the  repulsive  forces  are  not  Important  and  do  not 
significantly  deform  che  molecular  potentials.  This  Is  Che  case  for  asny  reac¬ 
tive  systems  that  exhibit  'harpooning*  behavior,  e.g. ,  alkali  and  metastsble  rare 
gas  collisions  with  halogens  and  Oj.  This  second  requirement  Is  a  mutter  of 
convenience,  since  the  general  effects  described  here  should  be  Important  even  In 
system  where  repulsive  forces  are  significant  In  the  region  of  electronic  state 
Interaction.  A  third  requirement  for  clear  observation  of  these  vibrational 
effects  is  the  need  for  the  collision  tlm  to  be  not  vastly  different  from  the 
vibrational  period  of  the  Intermediate  diatomic.  In  general  this  Is  the  case  for 
collision  energies  from  a  fraction  of  an  eV  to  several  hundred  eV,  depending  on 
the  syacam. 
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We  conclude  that  v  ^rational  timing  In  the  aubplcoaecoud  nag*  li  an  Important 
and  neglected  at  fact,  Che  aCudy  of  which  can  ravaal  the  lncarnal  aoclon  of  tran- 
alanC  acacaa  and  chui  provlda  datailad  Information  on  Cha  dynaaica  of  aolacular 
colllalona  and  cha  pocanclal  aurfacaa  involved.  ( 
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Ion-Beam  Neutralization-Reionization  Spectroscopy  of  Ion-Pair  Formation 
in  Reactions  of  He*(2JS)  and  He*(2'5)  with  02 

Thomas  M.  Miller1  and  Keith  T.  Gillen 
Molecular  Physics  Laboratory,  SRI  International,  Menlo  Park,  California  94025 
(Received  5  February  1980) 

For  He*  beams  neutralized  in  alkali  vapors  and  re  Ionized  in  Oj  gas,  information  on 
the  Intermediate  neutral  states  has  been  used  for  the  first  time  to  Isolate  the  differen¬ 
tial  cross  sections  for  ion-pair  formation  in  a  pair  of  reactions  occurring  simultane¬ 
ously.  The  reactions  of  He*(2%)  and  He*(2tS)  with  O,  are  strongly  coupled  by  a  shared 
Coulombic  product  channel,  and  evidence  exists  for  unique  three- surface  trajectories 
leading  to  ion-pair  products  He*  +  O," . 


At  ion  energies  E,  above  - 10  eV,  metastable 
rare -gas  beams1  can  be  produced2  efficiently  by 
charge-transfer  neutralization  of  the  parent  ion 
beam  in  an  alkali  vapor.  The  neutralized  beam 
(after  radiative  decay  of  short-lived  states)  is 
usually  an  unknown  mixture  of  the  two  metastable 
states  (2‘S  and  2*S  for  He*)  and  the  ground  state 
(*6^).  One  would  generally  expect  that  both  meas¬ 
urement  and  modification  of  the  neutral -beam 
composition  would  be  needed  for  isolating  the 
contributions  of  the  various  neutral-beam  com¬ 
ponents  to  subsequent  scattering  experiments. 
However,  we  have  used  only  a  knowledge  of  the 
states  produced  in  the  neutralization  reaction 

He  *  ♦  Rb  —  He*  ♦  Rb  *  (1) 

under  conditions  in  which  the  beam  composition 
is  not  determined,  to  resolve  and  investigate 
separately  the  ion-pair  reactions 

He*(2*S)  +Oa  -  He  *  +Oa'  (2) 

and 

He*(2‘SUOa-He*+Oa*.  (3) 

These  results  are  the  first  in  which  detailed  ion- 
pair -formation  data  have  been  obtained  for  two 
separate  entrance  channels  leading  to  a  common 
exit  channeL  We  will  show  that  Reactions  (2)  and 
(3)  have  remarkably  different  angular  distribu¬ 
tions  that  can  nevertheless  be  understood  from 
qualitative  estimates  of  the  deflection  functions 
for  the  important  trajectories  leading  to  ion-pair 
products.  Reaction  (3)  appears  to  be  the  first 
to  show  evidence  for  the  coupling  of  three  differ¬ 
ent  potential  surfaces  in  a  single  trajectory  lead¬ 
ing  to  ion-pair  products. 

Our  apparatus  and  basic  techniques  have  been 
described  in  detail.1,4  Briefly,  a  He*  beam  is 
produced  by  a  near -resonant  charge-transfer 
neutralization  of  a  fast  He  ’  beam  in  Rb  or  Na 
vapor.  Unreacted  ions  are  swept  aside  before 


the  beam  enters  a  collision  cell  filled  with  Oa. 

He  *  product  ions  from  this  second  collision  are 
detected  by  a  Channeltron  mounted  behind  an  127° 
electrostatic  energy  analyzer.  This  detector  sys¬ 
tem  can  be  rotated  around  the  scattering  center 
to  measure  the  product  He  *  angular  distributions. 
The  energy  analyzer  can  be  floated  at  various 
potentials  to  analyze  product  ions  with  high  or 
low  resolution.  A  separate  Channeltron  detector 
was  used  to  measure  the  He*  beam  energy  (£„) 
with  time-of-flight  (TQF)  techniques. 

For  a  He  *  beam  neutralized  to  He*(2*S)  in  Re¬ 
action  (1)  and  then  reionized  in  Reaction  (3),  the 
translational-energy  separation  EL  between  the 
original  He  *  beam  (£,)  and  the  final  He*  product 
(£')  is  given  bys 

Ex.  =  £1-£'=ff  +  /(Rb)-/(He*(21S)) 

+  /(He*(2,S))-A,(Oa)+£iol(0,'),  (4) 

where  K  is  the  Rb  contact  potential,  /[X]  is  the 
ionization  potential  of  X,  A,  is  tin  adiabatic  elec¬ 
tron  affinity,  and  E,tt< (Oa‘)  is  the  product  Oa* 
internal  energy.  The  first  three  terms  on  the 
right-hand  side  give  the  energy  shift  due  to  neu¬ 
tralization,  and  the  last  three  terms  are  for  the 
reionization  reaction.  The  third  and  fourth  terms 
cancel,  implying  that,  in  the  absence  of  radiative 
cascade,  any  other  intermediate  neutral  states 
yielding  the  same  Oa~  internal  states  in  the  re- 
lonlzation  step  would  produce  He4  product  ions 
of  Identical  energy- loss  values  Et.  In  contrast, 
neutralization  to  a  state  that  radiates  before  the 
reionization  collision  causes  the  third  and  fourth 
terms  in  Eq.  (4)  to  differ  by  the  energy  of  the 
photon(s)  radiated. 

A  recent  communication*  on  TOF  measurements 
of  excited-He-state  populations  produced  in  near- 
resr.iant  charge  transfer  of  He*  with  various 
alkali  vapors  demonstrates  that  100-200-eV 
charge  transfer  in  Rb  (or  Cs)  yields  significant 
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FIG.  1.  Representative  Ho*  product  energy  dlstrlbu- 
tloua  for  He*  beeme  formed  by  oharge  transfer  in  Rb 
(solid  line)  and  in  Na  (dashed  line). 


amounts  of  only  two  He  states:  2'S  and  2 3P.  The 
2sS  metastable  component  of  the  beam  comes  ex¬ 
clusively  from  rapid  radiative  decay  of  the  2 *P 
state.  Our  TOF  measurements  over  the  beam 
energy  range  30-200  eV  confirm  that  no  2’S  is 
produced  directly  in  Reaction  (1)  in  this  energy 
range.  Hence,  the  He  *  product  ions  from  Reac¬ 
tions  (2)  and  (3)  could  be  separated  by  1.14  eV 
(the  energy  of  the  2sp-2sS  photon)  in  accord  with 
Eq.  (4). 

Figure  1  presents  product  He*  laboratory  en¬ 
ergy  distributions  at  two  scattering  angles  for 
He*  collisions  at  He*  beam  energies  of 
-85.5  eV.  The  doubled-peaked  He  *  distributions 
result  from  He*  beams  produced  by  charge  trans¬ 
fer  in  Rb;  the  peak  separations  are  consistent 
with  the  predicted  shift  of  the  products  of  Reac¬ 
tion  (2)  to  lower  energy  due  to  the  2 3p  -  2 *S  pho¬ 
ton.  The  dashed  curves,  of  slightly  lower  ana¬ 
lyzer  resolution,  present  analogous  results  for 
a  beam  formed  in  Na  where  the  TOF  measure¬ 
ments  prove  that  He*(2,S)  is  the  only  major  beam 
component.  The  energy  widths  of  the  measured 
He*  profile  for  both  Reactions  (2)  and  (3)  are  con¬ 
volutions  of  the  beam  width  (typically  <0.5  eV 
full  width  at  half  maximum),  the  detector  band¬ 
pass  (adjustable,  but  typically  0.5  eV  at  highest 
resolution  used),  and  the  widths  of  the  product 
state  distributions  at  the  measured  angle.  Clear¬ 
ly  these  product  state  distributions  are  narrow 
enough  that  the  two  reactions  can  be  easily  re¬ 
solved. 


6  (degrees) 


FIG.  2.  Product  He*  laboratory  angular  distributions 
from  reactions  of  He*(21S)  and  He*(2%)  with  O, . 


Similar  energy  scans  at  other  angles  demon¬ 
strate  that  these  peaks  neither  broaden  nor  shift 
in  energy  over  the  angular  range  in  which  there 
is  significant  scattered  intensity.  This  result 
suggests  a  particularly  simple  approach  for  meas¬ 
uring  the  angular  distributions  of  the  products  of 
Reactions  (2)  and  (3)  quickly  and  separately. 

After  He*  production  with  Rb  charge  transfer, 
separate  angular  distributions  are  determined 
with  the  energy  analyzer  at  high  resolution  and 
set  to  transmit  ions  with  energies  near  the  peak 
of  the  signal  for  Reaction  (2)  and  for  Reaction  (3). 
Reaction  (2)  can  be  studied  independently  for  a 
He  *(2*5)  beam  produced  in  Na,  and  energy-inte¬ 
grated  angular  distributions  are  found  to  be  in 
excellent  agreement  with  those  obtained  with  the 
simpler  method.  The  product  He*  angular  dis¬ 
tributions  shown  in  Fig.  2  are  typical  results  that 
demonstrate  dramatic  differences  between  the 
two  reactions. 

Analysis  of  several  recent  measurements  of 
detailed  differential  cross  sections  for  ion-pair 
formation  has  yielded  much  insight  into  the  reac¬ 
tion  dynamics.7  The  projectile  is  usually,  but  not 
always,4  a  ground-state  alkali  atom.  All  these 
systems4,7  are  strongly  Influenced  by  crossings 
at  large  internuclear  distances  (typically  3-7  A) 
between  an  essentially  flat  Incoming  covalent  sur¬ 
face  and  a  strongly  attractive  Coulombic  surface 
asymptotically  a  few  electronvolts  higher.  Tran¬ 
sitions  can  occur  as  this  crossing  region  is  tra¬ 
versed.  At  collision  energies  above  a  few  elec¬ 
tronvolts,  the  products  can  separate  as  ion  pairs 
with  a  very  large  reaction  cross  section. 

In  collisions  involving  metastable  rare-gas 
atoms,4  the  interacting  surfaces  are  often  also 


F-2 


777 


Volumi  44,  NuMin  12 


PHYSICAL  REVIEW  LETTERS 


24  Makch  1980 


8  (A) 


FIG.  3.  Schematic  He*-0,  potentials,  showing  all 
the  possible  surface  symmetries  In  the  lowest  sym¬ 
metry  (C, )  orientation.  All  three  asymptotic  states 
are  coupled  by  3A"  surfaces. 


Imbedded  with  a  continuum,  and  competition 
from  Penning  ionization  and  excitation  transfer 
reactions  may  influence  the  ion-pair  distributions. 
An  additional  and  unique  feature  of  Reactions  (2) 
and  (3)  is  their  expected  direct  coupling  to  each 
other,  as  indicated  in  Fig.  3,  which  shows  the 
interaction  of  He*  and  O,  at  a  fixed  0-0  distance 
(nominally  the  O,  equilibrium  bond  length).  The 
symmetries  of  the  various  surfaces  in  C,  geom¬ 
etry  are  indicated  for  both  incoming  reaction 
pairs  and  for  the  ionic  product  channel.  Note 
that  both  reactants  can  couple  to  the  ion-pair 
channel;  note  especially  that  the  3A  '  surfaces 
provide  a  mechanism  for  coupling  of  all  three 
asymptotic  states. 

The  angular  position  of  the  three  distinct  peaks 
(at«0°,  2.6°,  and  £5°)  evident  in  Fig.  2  for  the 
2*s  reaction  and  the  two  peaks  (near  2s  and  4.2s) 
for  reaction  of  He*(2*S)  can  now  be  understood 
qualitatively  from  an  examination  of  the  possible 
trajectories  in  Fig.  3  that  yield  ion -pair  products. 
For  the  He*(2*S)  +0,  reactants,  the  two  paths 
“aedc"  (called  the  “ionic"  path)7  and  “afdc”  (the 
“covalent"  path)  are  possible.  Two  prominent 
peaks  would  be  expected  in  the  angular  distribu¬ 
tion.  The  peak  at  larger  angles  (from  more 
attractive  Interactions)  would  be  a  rainbow  from 
the  portion  ot  the  deflection  function  associated 
with  trajectories  following  the  ionic  path.  A  sec¬ 
ond  maximum  at  smaller  angles  is  due  to  the 
attractive  Interactions  on  the  outgoing  portion  of 
the  covalent  trajectories.  This  expectation  is 
based  on  analogous  results  for  many  other  ion- 
palr  reactions.4'7 


The  He*(2*S)+0*  system,  however,  is  the  first 
one  studied  in  which  the  reactant  channel  is  not 
the  lowest  coupled  state  asymptotically.  In  addi¬ 
tion  to  the  “ionic"  path  ( bdedc )  and  the  “covalent” 
path  ( bgc )  to  the  ion-pair  products,  there  is  a 
three-surface  trajectory  (bdfdc).  The  rainbow 
peak  from  the  ionic  trajectory  would  be  expected 
to  be  at  a  larger  angle  than  the  corresponding 
rainbow  for  the  2 3S  reaction,  since  more  attrac¬ 
tion  occurs  during  the  trajectory  (an  extra  tra¬ 
versal  of  region  d).  Likewise,  the  covalent  peak 
from  Reaction  (3)  would  be  expected  at  smaller 
angles  than  that  from  Reaction  (2).  The  three- 
surface  trajectory  could  then  be  responsible  for 
the  prominence  of  the  intermediate  peak  at  -2.6°. 
Qualitative  estimates  of  the  deflection  functions 
expected  for  the  five  different  reaction  paths  are 
then  in  good  agreement  with  the  measured  rela¬ 
tive  angular  positions. 

Better  estimates  of  the  product  distributions 
require  deflection  function  calculations  that  at¬ 
tempt  to  match  the  data  quantitatively  by  account¬ 
ing  for  surface  coupling  parameters  and  branch¬ 
ing  ratios  and  effects  from  competition  with  other 
channels  (including  continuum  coupling).  Classi¬ 
cal  trajectory  calculations  confirm  and  extend 
the  general  conclusions  presented  here;  in  par¬ 
ticular,  the  three-surface  peak  at  2.6°  Is  enhanced 
by  trajectories  (bdc)  that  do  not  nominally  reach 
the  second  crossing  and  is  modified  further  by 
strong  deviations  from  diabetic  behavior  at  the 
inner  crossing.  These  calculations  and  the  de¬ 
tailed  experimental  results  will  be  presented 
elsewhere.' 
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